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L  INTRODUCTION 


Modem  naval  warfare,  especially  undersea  warfare,  depends  heavily 
on  the  vessel  being  acoustically  silent.  A  major  source  of  radiated  noise  is 
the  vibration  of  shipboard  components.  The  reduction  of  these  vnbrations  is 
of  utmost  importance  if  a  ship  is  to  accomplish  its  mission.  One  method  of 
vibration  damping  that  shows  promise  in  damping  over  a  broad  spectrum 
of  low  frequency  vibration  is  constrained  viscoelastic  layer  damping.  The 
constrained  viscoelastic  layer  method  uses  the  high  energy  dissipation 
characteristic  of  viscoelastic  materials  during  periodic  motion  of  shear 
deformation  to  absorb  and  dissipate  the  vibrational  energy  of  the  system  in 
question.  Unfortunately  the  design  and  analysis  of  such  constrained 
viscoelastic  layer  systems  is  difficult,  due  in  part  to  the  following: 

•  The  material  properties  of  viscoelastic  damping  materials  vary 
greatly  with  temperature  and  frequency. 

•  Closed  form  solutions  to  the  equations  of  motion  for  constrained 
layer  system  exist  only  for  beams  and  plates  with  simple 
boundary  conditions. 

•  The  exact  complex  valued  eigenvalue  analysis  for  constrained 
viscoelastic  layer  damping  systems  using  the  finite  element 
method  requires  large  amounts  of  computer  storage  and  CPU  time. 
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Johnson  and  Kienholz  developed  the  Modal  Strain  Energy  (MSE) 
method  which  uses  the  ratio  of  strain  energy  for  each  mode  shape  to 
approximate  the  modal  damping  of  a  structure  for  a  given  constrained 
viscoelastic  damping  system  [Ref.  1].  This  method  is  very  attractive  due  to 
its  simple  concept  and  very  useful  because  it  can  be  applied  to  any  general 
cases  with  arbitrary  shape  by  using  the  finite  element  method.  However, 
its  effectiveness  compared  with  experiments  were  reported  for  only  a  few 
cases.  Maurer  examined  the  effectiveness  of  the  MSE  method  for  two 
damped  plate  configurations;  1)  a  simple  sandwich  configuration,  and  2) 
a  plate  with  a  milled  pocket  with  damping  material  inserted  and  a  welded 
cover  plate  acting  as  a  constraining  layer  in  a  previous  work  [Ref.2]. 
However,  he  could  not  verify  for  the  milled  pocket  plate  case  since  the  cover 
plate  warped  and  delaminated  itself  from  the  damping  material  during 
welding,  resulting  in  negligible  damping  [Ref.  2]. 

For  certain  naval  applications  the  components  to  be  damped  will  be 
thick  in  construction  and  may  be  exposed  to  an  unfriendly  environment. 
Therefore,  in  this  research  the  pocket  plate  configuration  and  a  second 
milled  plate  using  a  “floating  element”  in  conjunction  with  constrained 
layer  damping  are  investigated.  These  damping  treatments  are  compared 
with  the  simple  sandwich  type  treatments  consisting  of  single  and  double 
constraining  layers.  Simple  plate  geometries  were  used  to  facilitate  the 
experimental  and  computational  effort.  Therefore,  this  paper  addresses 
the  experimental  testing  and  analysis  of  four  thick  aluminum  plates,  each 
with  a  different  constrained  viscoelastic  layer  damping  treatment.  The 
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effectiveness  of  the  MSE  method  is  investigated  for  each  different  damping 
treatment  by  evaluating  its  accuracy  in  modal  damping  value  prediction 
compared  with  experimental  results  and  usefulness  as  a  possible  design 
tool. 
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II.  THEORY 


A.  VISCOELASTIC  MATERIAL 

Viscoelastic  materials  of  interest  for  general  naval 
applications  are  polymeric  compounds  made  up  of  long  molecular  chains. 
These  molecular  chains  can  be  strongly,  or  weakly,  linked  together, 
depending  on  their  chemical  composition  and  processing.  The  damping 
characteristics  of  viscoelastic  arise  from  the  deformation  and  recovery  of 
the  polymer  network.  Material  properties  of  a  viscoelastic  material  vary 
with  temperature  and  frequency.  As  such,  the  damping  characteristics  of 
a  system  will  vary  as  its  operating  environment  changes.  [Ref  3] 

Temperature  will  have  the  greatest  effect  on  the  material  properties 
of  damping  materials  [Ref  3].  This  effect  is  shown  in  Figure  2.1,  where 
four  distinct  regions  are  observed.  The  lowest  temperature  region  is  the 
glassy  region  where  the  material’s  storage  modulus  is  at  its  maximum 
value,  and  the  loss  factor  is  at  a  minimum.  In  the  glassy  region  the 
modulus  decreases  slowly  with  temperature  increase,  whereas  the  loss 
factor  increases  rapidly  with  temperature.  The  second  region  is  the 
transition  region  where  the  modulus  decreases  rapidly  with  increasing 
temperature  and  the  loss  factor  reaches  its  peak  value.  The  third  region  is 
the  rubbery  region  where  both  the  modtilus  and  loss  factor  are  at  low  values 
and  show  little  variation  with  temperature.  The  fourth,  and  last,  region  is 
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the  flow  region  and  characterizes  the  behavior  of  some  materials,  mostly 
ceramics,  at  high  temperatures.  It  should  be  noted  that  the  transition 
region  may  vary  in  width  from  20  up  to  a  width  of  200  °C.  [Ref  3] 

The  effect  of  frequency  on  viscoelastic  materials  is  not  as  great  as  that 
of  temperature.  The  modulus  of  the  viscoelastic  always  increases  with 
increasing  frequency.  The  loss  factor  will  initially  increase  with  frequency, 
then  peak,  and  subsequently  decrease  as  frequency  increases.  A  plot  of 
storage  modulus  and  loss  factor  versus  frequency  is  shown  in  Figure  2.2.  It 
should  be  noted  that  this  plot  is  over  a  range  of  approximately  ten  decades, 
and  hence  it  becomes  obvious  that  a  temperature  change  of  a  couple  degrees 
will  have  a  much  greater  effect  on  damping  than  a  minor  change  in 
frequency.  [Ref  3] 

Linear  viscoelastic  materials  behave  in  a  hysteretic  manner  under 
cyclic  excitation.  Therefore,  the  mechanical  behavior  of  a  viscoelastic 
material  during  steady  state  vibration  is  best  described  by  using  a  complex 
stiffness,  k*  [Ref  3]. 

k*  =  k(l  +  iq)  (2.1) 

where, 

q  =  material  loss  factor 

The  use  of  a  complex  stiffness  then  leads  to  the  use  of  a  complex  Yoimg’s 
modulus  and  shear  modulus  [Ref  3]. 
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E*  =  E(1  +  iri) 
G*  =  G(1  +iri) 


(2.2) 

(2.3) 


This  concept  of  the  complex  modulus  is  used  in  subsequent  analysis. 

Viscoelastic  material  properties  are  commonly  displayed  using  a 
“reduced  frequency  nomogram.”  The  reduced  frequency  nomogram 
displays  the  variation  of  the  viscoelastic  material’s  loss  factor  and  modulus 
with  temperature  and  frequency.  The  “reduced  frequency”,  fat.  Is  an 

empirically  determined  function  that  accounts  for  the  viscoelastic’s 
temperature  and  frequency  dependence,  and  allows  data  for  wide  range  of 
temperature  and  frequencies  to  be  plotted  on  the  same  graph  [Ref  4].  The 
reduced  frequency  nomogram  for  3M  ISD  -  112  is  shown  in  Figure  2.3.  To 
find  the  loss  factor  and  modulus  using  the  nomogram,  enter  with  the 
desired  temperature  and  frequency.  Follow  the  frequency  line  horizontally 
and  the  temperature  line  diagonally  down  the  page  until  the  two  intersect. 
Then  go  vertically  up  or  down  to  intersect  the  shear  modulus  or  loss  factor 
curves.  Finally,  read  the  value  of  the  shear  modulus  or  loss  factor 
horizontally  from  the  scale  on  the  left  [Ref.  4]. 


R  CONSTRAINED  VISCOELASTIC  LAYER  DAMPING 

A  simple  constrained  layer  damping  treatment  consists  of  a  base 
layer  (the  structure  to  be  damped),  a  damping  layer,  and  the  constraining 
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layer.  This  configuration  is  shown  in  Figure  2.4  with  the  thicknesses  of  the 
damping  and  constraining  layers  exaggerated  for  clarity. 

The  physical  mechanism  of  damping  can  be  explained  by  referring  to 
Figure  2.4.  When  the  base  layer  is  deformed  in  a  mode  of  vibration,  the 
surface  away  from  the  neutral  axis  elongates,  stretching  the  viscoelastic 
material.  The  top  layer,  being  a  stiff  elastic  material,  tends  not  to  elongate, 
and  thereby  “constrains”  the  viscoelastic  material.  Consequently,  the  cyclic 
motions  of  vibration  induce  a  cyclic  shearing  strain  in  the  viscoelastic. 
This  cyclic  shearing  strain,  together  with  its  associated  hysteresis  loop 
cause  the  vibrational  energy  to  be  dissipated  as  heat.  For  the  constraining 
layer  to  be  effective,  its  stiffness  should  not  exceed  that  of  the  base  layer, 
[Ref.  4  &  5] 


C.  SYSTEM  EQUATIONS  OF  MOTION 

Continuous  systems,  such  as  plates,  possess  distributed 
characteristics  of  mass,  damping,  and  stiffness.  Classical  vibration 
analysis  of  such  systems  involves  the  formation  of  a  mathematical  model 
that  discretizes  the  system  into  a  finite  number  of  components  in  order  to 
approximate  the  total  system.  Such  a  formulation  results  in  the  following 
equation: 

[M]  {x(t))  +  [C]  {x(t))  +  [K]  {x(t)}  =  {F(t)}  (2.4) 
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where, 


[M]  =  system  mass  matrix 
[C]  =  system  damping  matrix 
[K]  =  system  stiffness  matrix 
{F(t)}  =  external  excitation  vector 
{x(t)}  =  displacement  vector 

For  an  undamped  system  without  excitation,  the  above  equation  reduces  to 
the  eigenvalue  problem. 

;m' +  ;K;!x(t)i  =  0  (2.5) 

This  equation  is  then  transformed  to  modal  space  using  the  linear 
transformation: 

{x(t)}  =[<p]  {q(t)}  (2.6) 


where, 

[<{)]  =  modal  matrix 

{q(t)}  =  modal  response  vector 

Using  this  linear  transformation  the  equation  of  motion  can  then  be  solved 
for  the  undamped  modal  frequencies  and  mode  shapes. 

To  solve  for  the  frequency  response  of  a  damped  system,  the  linear 
transformation  is  applied  to  equation  (2.4): 

[M]  M  {q(t)}  +  [C]  W  {q(t))  +  [K]  M  {q(t)}  =  {F(t)}  (2.7) 
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Assuming  that  the  damping  matrix  [C]  is  proportional  to  a  linear 
combination  of  the  stiffness  matrix  [K]  and  mass  matrix  [M],  the  damping 
matrix  can  then  be  diagonalized  using  the  same  linear  transformation 
used  to  diagonalize  [K]  and  [M]  in  equation  (2.5)  above.  The  diagonal  terms 
of  the  damping  matrix  then  become  (TliCOj),  where  qj  equals  the  modal  loss 
factor  and  coi  is  the  natural  frequency  of  the  i^^  mode  [Ref.  1]  .  Using  this 
approximate  diagonal  damping  matrix  results  in  a  system  of  uncoupled 
modal  equations  of  motion: 

qj(t)  +  q,cq4(t)  +  ofqj(t)  =  fi(t)  (2.8) 

where, 

qi(t)  =  modal  acceleration  of  i^^  mode 
qi(t)  =  modal  velocity  of  i^^  mode 
q(t)  =  modal  displacement  of  i^^  mode 
Hi  =  modal  loss  factor  of  i^^  mode 
cq  =  i^^  natural  frequency 
fi(t)  =  modal  force  in  mode 

j  =  VT 

Assuming  that  a  sinusoidal  excitation  produces  a  sinusoidal  response, 

{fi:t)}  =  (flei‘^  {q(t)}  =  {Q)ej“‘  (2.9) 

the  response  for  the  i^^  mode  is  then  solved  to  be: 
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(2.10) 


-  a>2  +  jcojriiCi) 

Subsequently,  the  response  of  the  physical  system  can  be  found  using: 

(x(t)}=[(!)]{Q}6j‘^  (2,11) 

D.  MODAL  STRAIN  ENERGY  METHOD 

The  equations  of  motion  used  to  define  the  response  of  a  system  with 
viscoelastic  materials  need  a  complex  eigenvalue  analysis.  However,  the 
actual  solution  of  these  equations  may  be  quite  difficult.  This  is  especially 
true  when  the  system  to  be  analyzed  is  comprised  of  materials  whose 
properties  vary  with  both  temperature  and  frequency.  Finite  element 
techniques  are  generally  used  to  compute  the  response  of  complicated 
systems.  However,  for  the  case  of  varying  material  properties  many  time 
consuming  and  costly  runs  must  be  made  with  the  material  properties 
changing  at  each  frequency  increment  [Ref.  1].  In  addition  to  the  costly 
analysis  of  a  single  design  configuration,  changes  in  design  options,  design 
requirements,  of  the  search  for  an  optimum  design  can  make  the  expense 
of  finite  element  analysis  too  great.  The  development  of  the  Modal  Strain 
Energy  (MSE)  method  by  Johnson  and  Kienholz,  however,  makes  the  finite 
element  analysis  of  complex  viscoelastically  damped  structures  a  viable 
option  [Ref.  1]. 
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The  MSE  method  assumes  tiiat  a  damped  structure  can  be 
represented  by  the  normal  modes  of  the  associated  undamped  system  if 
appropriate  damping  terms  are  inserted  into  Uie  uncoupled  modal 
equations  of  motion  (eqn  2.8)  [Ref.  1].  The  MSE  method  further  assumes 
that  the  damping  matrix  of  equation  (2.4)  can  be  diagonalized  by  the  same 
real  modal  matrix  that  diagonalizes  the  system  mass  and  stiffness 
matrices  of  equation  (2.4). 


The  modal  loss  factors  of  equation  (2.8)  are  calculated  using  the 
undamped  mode  shapes  and  material  loss  factors  for  each  material.  Since 
the  material  loss  factors  of  the  structure  to  be  damped  and  the  constraining 
layers  are  quite  small  compared  to  those  of  common  viscoelastic  cores,  the 
modal  loss  factors  can  be  estimated  using  the  following  equation  [Ref.  1]: 


=  h 


(dY^ 


(2.12) 


whore, 

=  system  loss  factor  in  the  r^^  mode 

Ti^v^  =  material  loss  factor  of  viscoelastic  core  at  the 
natural  frequency 

Vy  ^  *=  strain  energy  in  the  viscoelastic  core  at  r^^  mode 
t:  strain  energy  in  the  entire  structure  at  r^^  mode 
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The  modal  strain  energies  can  be  obtained  for  finite  element  analysis,  and 
are  a  standard  output  option  of  the  NASTRAN  finite  element  code  [Ref  7]. 

The  modal  frequency  response  of  the  structure  is  then  calculated 
using  the  modal  loss  factors  found  in  equation  (2.12).  When  computing  the 
modal  frequency  response  of  a  damped  structure,  the  modal  properties  in 
the  system  matrices  are  assumed  to  be  constant.  However,  viscoelastic 
materials  have  storage  moduli  which  are  frequency  dependent.  To  account 
for  this  frequency  dependence,  Johnson  and  Kienholz  devised  the  following 
correction  factor  to  be  applied  to  the  modal  loss  factors  calculated  in 
equation  (2.12)  [Ref  8]. 


(2.13) 


where, 

=  corrected  modal  loss  factor  at  the  r^^  mode 
G2(fr)  =  viscoelastic  shear  modulus  at  the  r*-^  modal 
frequency 

(^,ref  =  reference  viscoelastic  shear  modulus  used  in  the 
frequency  response  calculation 
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Figure  2.1.  Variation  of  Viscoelastic  Material 

Properties  with  Temperature  [Ref.  3]. 


Figure  2.2.  Variation  of  Viscoelastic  Material 
Properties  with  Frequency  [Ref.  3]. 
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Figure  2.3.  Temperature  Frequency  Nomogram  for  3M  ISD  -  112. 
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Figure  2.4.  Single  Constrained  Layer  Configuration. 


IIL  DESIGN  OF  DAMPED  PLATES 


A .  GENERAL  SPECIMEN  CONFIGURATIONS 

For  experimental  testing  and  analysis  purposes,  four  different,  yet 
related,  constrained  layer  damping  treatments  were  selected  in  addition  to 
an  undamped  “reference”  plate.  Two  of  the  damping  configurations  were 
simple  sandwich  treatments  consisting  of  one  and  two  viscoelastic  layers 
respectively.  Another  treatment  was  the  “pocket  plate”  which  was 
previously  investigated  by  Maurer  [Ref  2].  The  pocket  plate  was  made  from 
a  solid  plate  which  was  then  milled  to  accept  damping  material  and  a  cover 
plate.  The  final  damping  treatment  was  a  “floating  element”  configuration 
consisting  of  a  solid  plate  milled  to  accept  a  double  layer  of  damping 
material  and  a  welded  cover  plate.  Section  views  of  these  damping 
configurations  are  shown  in  Figure  3,1. 

The  purpose  of  the  pocket  plate  is  to  protect  the  viscoelastic  material 
from  materials  such  as  oil  and  salt  water,  which  may  harm  the 
viscoelastic.  Since  previous  attempts  at  using  a  welded  cover  plate  were 
unsuccessful  due  to  the  heat  of  welding  causing  a  delamination  between 
the  damping  material  and  cover  plate  [Ref.  2],  it  was  decided  to  move  the 
viscoelastic  material  from  the  welding  point  by  recessing  it  into  a  shallow 
pocket  of  its  own  as  shown  in  Figure  3.1. 


The  "floating  element”  concept  evolved  because  the  welded  cover  plate 
of  the  pocket  plate  configuration  does  not  produce  the  damping  reaction  that 
a  true  constraining  layer  would  provide.  If  the  cover  plate  is  welded  to  the 
surrounding  structure  it  cannot  deform  in  bending  as  much  as  an 
iinwelded  constraining  layer,  thereby  causing  a  reduction  in  the  damping 
capability  of  the  system.  By  using  a  piece  of  metal  in  the  milled  pocket  with 
dimensions  slightly  smaller  than  the  surrounding  pocket,  along  with  two 
layers  of  viscoelastic  and  a  welded  cover  plate,  a  true  constraining  layer 
effect  should  be  obtained.  This  configuration  is  shown  in  Figure  3.1. 

In  order  to  approximate  a  possible  system  to  be  damped,  a  plate 
with  large  dimensions  was  selected.  The  dimensions  of  the  plates  used  for 
the  damping  treatments  are  114.3  cm  (45  in)  in  length  and  38.1  cm  (15  in) 
in  width.  In  addition,  possible  naval  applications  for  this  type  of  damping 
would  probably  consist  of  thick  plate  members.  For  this  reason  it  was  also 
decided  that  thick  plates  would  be  used  for  the  damping  treatments.  All  the 
base  plates  and  constraining  layers  were  made  of  a  standard  6061— T6 
aluminum  alloy.  The  design  and  selection  of  viscoelastic  layer,  base  layer, 
and  constraining  layer  thicknesses  is  discussed  in  the  following  sections. 


R  DESIGN  OF  THE  SINGLE  DAMPING  LAYER  CONFIGURATION 

In  an  attempt  to  approximate  system  loss  factors  and  hence 
determine  viscoelastic  and  constraining  layer  thicknesses  for  maximum 


damping,  a  method  developed  by  Nashif  [Ref.  9]  based  on  an  analysis  of 
simple  sandwich  damping  systems  by  Ross,  Kerwin,  and  Ungar  [Ref.  10] 
was  usea.  The  Ross  -  Kerwin  -  Ungar  (RKU)  equations  are  base  of  the 
analysis  of  the  simple  sandwich  system  shown  in  Figure  3.2. 

To  find  the  loss  factors  of  the  damped  system,  the  flexural  rigidity  of 
the  system  must  first  be  determined.  For  the  above  system,  the  flexural 
rigidity,  El,  is  written  as  [Ref.  9]: 


EI  = 


E,H^ 


1  ,  E2H^ 
12 


+ 


EsH|  + 


E,H,D‘ 


12  12  12 
+  E2H2(H2i  -  D)2+  EaHalHai  -  Df 

-|  ^  +  ^|^H2i  -  D)  +  EaHsiHsi  -  D) 


•  D 
1+g 


(3.1) 


where, 

+  g(E2H2H2l  +  E3H3H31) 

F  R  -  ^3.2) 

EiHi  +  +  g(EiHi  +  E2H2  +  E3H3) 


(3.3) 

(3.4) 

g _ 

E3H3H2K2 

(3.5) 

E  =  Young’s  modulus 
G  =  Shear  modulus 


I  =  second  moment  of  area 
H  =  thickness  of  member 
K  =  wave  number 


Subscripts  refer  to  the  layers  labeled  in  Figure  3.1.  No  subscript  refers  to 
the  composite  system 


For  a  simply  supported  plate  the  wave  numbers  and  modal 
frequencies  are  found  using  [Ref.9]: 


—  iMim' 


EH^ge 


12(l-v2)Hp 


(3.6) 


(3.7) 


where, 

a  =  semi-wave  length  of  the  plate 
b  =  semi-wave  width  of  the  plate 
V  =  Poisson’s  ratio  of  the  composite  plate 
p  =  density  of  the  composite  plate 
gc  =  gravitational  constant 

To  introduce  damping  into  the  equations  it  is  necessary  to  use  the 
complex  modulus  concept  expressed  in  Section  11.  Substituting  the 

appropriate  expressions  for  the  complex  shear  and  Young’s  modulus  into 
equations  (3.1),  and  assuming  that  damping  in  the  base  layer,  (tij),  is 

small,  and  that  the  extensional  stiffness  of  the  damping  layer  is  small 
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(since  E2  «  and  E2  «  E3),  the  following  expressions  can  be  arrived  at 
[Ref.9]: 


EH^  =  EiH?  +  E3H?  +  -1^-  (a  -  U  -  dk  (3.8) 

c2  +  d^ 

EH^-n  =  EaHana  +  -  ^  (a  -  (3  -  3)im  (3.9) 

c2  +  d^ 


where, 

a  =  gEiHiE3H3H§i{c(l  -  •n2n3)  +  d(ri2  +  ns)  +  j[c(Ti2  +  Ha)]}  (3.10) 

p  =  E1H1E2H2H31  [c  +  dTi2  +  j(cTi2  -  d)]  (3.11) 


a 


=  2gE2H2E3H3H2iH3i 


I  c(l-2n2Tl3-r|i)  +  d(2Tl2  +  Tl3-ri|Tl3)  \ 

1+  j[c(2Ti2  +  n3  -  nina)  -  d(l  -  2Ti2ri3  -  Tli)]/ 


(3.12) 


c  =  EiHi(l  +  g)  +  gE3H3{l-ri2Tl3)  (3.13) 

d  =  gEiHiri2  +  gE3H3(Ti2  +  ns)  (3.14) 

j  =  iTT  (3.15) 


These  equations  were  then  applied  to  estimate  the  loss  factors  of  the  simple 
three— layer  sandwich  plate.  The  equations  can  be  simplified  by  assuming 
that  damping  in  the  constraining  layer  (T13)  is  negligible  [Ref.9]. 


Since  the  boundary  conditions  for  the  plate  used  in  this  research 
(free-free-free-free)  do  not  correspond  to  the  simply  supported  conditions  on 
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which  equation  (3.7)  is  based,  modal  frequencies  for  the  free-free  case  were 
estimated  using  results  from  finite  element  analysis  of  a  free-free  plate  and 
equation  (3.6).  The  natural  frequencies  of  an  undamped  plate  with  the 
dimensions  previously  given,  and  a  thickness  of  1.91  cm  (0.75  in)  were 
found  using  a  normal  mode  extraction  in  NASTRAN.  By  substituting  these 
modal  frequencies  into  equation  (3.6)  and  estimate  of  the  wave  parameter 
for  each  mode,  K^,  was  obtained.  Then,  by  using  an  iterative  procedure 
outlined  in  Reference  [9],  modal  loss  factors  were  estimated  for  different 
layer  thicknesses  over  a  temperature  range  of  0.0  °C  to  37.8  °C  (30  to 
100  °F). 


The  previous  equations  are  easily  programmed  to  compute  loss 
factors  for  a  wide  variety  of  conditions.  The  variation  of  viscoelastic 
material  properties  with  temperature  and  frequency  was  accounted  for 
using  a  curve-fit  to  the  reduced  frequency  nomogram  developed  by  Drake 
[Ref.  11].  The  material  data  for  the  following  curve-fit  equations  is  from  the 
University  of  Dayton  Research  Institute  [Ref.  12]. 


logio(M)  =  logio(ML) 


2  logioj^ 


ROM 

ML 


1  .  (FROMF 
FR 


(3.16) 


logio(ETA)  =  logiolETAFROL) 

^[a(sl  +  sh)  +  {sl-  sni  1  -  h  +  A^)] 

(3.18) 
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(3.19) 


^  _  logio(FR)  -  logio(FROL) 

A=  ^ 

where, 

M  =  viscoelastic  modulus 
ETA  =  viscoelastic  material  loss  factor 
FR  =  reduced  frequency  (Hz) 

F  =  frequency  (Hz) 

T  =  temperature  (°F) 

and, 

T0  =40  °C(104  °F) 

FROM  =  2.0x104  Hz 
MROM  =  4.75x106  Pa  (688.94  psi) 
n  =  0.275 

ML  =  6.0x104  Pa  (8.7  psi) 

ETAFROL  =  1.08 
SL  =  0.45 
SH  =  -0.55 
FROL  =  5000  Hz 
C  =  2.5 

In  addition  to  the  above  constants,  a  Poisson’s  ratio  of  0.49  and  a  density  of 
0.909  gram  per  cubic  centimeter  (0.035  Ibm/in^)  was  used  for  ISD-112 
[Ref.  12].  The  following  material  properties  were  used  for  6061-T6 
aluminum  [Ref.  13]. 
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E  =  70  GPadOxlO^  psi) 

V  =  0.33 

p  =  2.7  ^m/cm^  (0.0968  Ibm/in^) 

Using  the  previous  equations  and  material  properties,  a  computer 
program  was  written  to  compute  estimated  modal  loss  factors  and 
frequencies  for  a  variety  of  base  layer,  viscoelastic  layer,  and  constraining 
layer  thicknesses.  A  listing  of  this  program  appears  in  Appendix  A. 
Modal  loss  factors  were  computed  for  base  layer  thicknesses  of  9.53  mm 
(0.375  in)  to  19.05  mm  (0.75  in)  in  3.18  mm  (0.125  in)  increments.  For  each 
base  layer  thickness,  the  viscoelastic  thickness  was  varied  from  0.38  mm 
(0.015  in)  to  1.52  mm  (0.060  in)  in  0.38  mm  increments,  and  the 
constraining  layer  thickness  was  varied  from  1.59  mm  (0.0625  in)  to  6.35 
mm  (0.25  in)  in  1.59  mm  increments.  In  addition,  loss  factors  were  also 
computed  for  a  viscoelastic  thickness  of  0.127  mm  (0.005  in).  From  the 
results  of  the  analysis,  a  carpet  plot  [Ref.  3],  was  made  for  each  of  the  base 
layer  conditions.  The  carpet  plot  reflects,  for  the  the  first  mode,  the 
maximum  loss  factor  and  its  corresponding  temperature  for  each 
viscoelastic  layer/constraining  layer  thickness  configuration.  The  carpet 
plot  for  a  base  layer  thickness  of  12.7  mm  (0.50  in)  is  shown  in  Figure  3.3. 
Based  on  the  carpet  plots  and  a  desire  for  maximum  damping,  as  well  as  a 
system  which  could  be  moved  easily,  a  base  layer  thickness  of  12.7  mm(0.50 
in),  a  viscoelastic  thickness  of  0.38  mm  (0.015  in),  and  a  constraining  layer 
thickness  of  6.35  mm  (0.25  in)  was  selected.  The  total  system  thickness  was 
approximately  19.05  mm  (0.75  in).  This  total  system  thickness  would  be 
maintained  for  all  subsequent  damping  configurations. 
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To  maintain  continuity  between  damping  systems,  the  milled  “pocket 
plate”  was  given  a  viscoelastic  thickness  of  0.38  mm  and  a  cover  plate 
thickness  of  6.35  mm  for  a  total  system  thickness  of  19.05  mm  (0.75  in).  In 
order  to  keep  the  heat  of  welding  away  from  the  viscoelastic  material,  the 
ISD— 112  was  recessed  into  a  shallow  pocket  as  indicated  previously  in 
Figure  3.1,  and  as  shown  in  the  pocket  plate  system  arrangement  in 
Figure  3.4.  Detail  drawings  of  the  pocket  plate  are  shown  in  Appendix  C. 


C.  DESIGN  OF  THE  DOUBLE  DAMPING  LAYER 

The  design  of  the  double  layer  damping  system  was  accomplished 
using  a  modification  of  the  RKU  analysis  used  in  the  previous  section.  The 
RKU  equations  are  used  by  working  from  the  top  layer  of  the  damping 
system  down  towards  the  base  layer.  As  shown  in  Figure  3.5,  the  H3' 
constraining  layer  along  with  the  H2'  viscoelastic  layer  are  combined  with 
the  HI'  layer  to  form  a  three-layer  system.  Using  the  RKU  equations,  the 
stiffness  of  this  system  is  computed  and  considered  to  be  the  equivalent 
stiffness  of  the  top  three  layers  of  the  total  constrained  layer  damping 
system.  The  top  three  layers  were  then  considered  as  a  single  layer  with 
the  equivalent  stiffness  previously  calculated,  and  the  RKU  equations  were 
again  applied  to  compute  estimated  modal  loss  factors  for  the  entire  double 
layer  damping  system  [Ref  3]. 
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To  maintain  continuity  among  all  the  damping  configurations,  a 
base  layer  thickness  of  12.77  mm  (0.50  in)  was  chosen,  and  a  total  system 
thickness  of  19.05  mm  (0.75  in)  was  maintained.  A  design  for  high 
damping  was  then  selected  by  computing  modal  loss  factors  for  the 
constraining  and  damping  layer  thickness  combinations  shown  in  Table 
3.1.  The  estimated  modal  loss  factors  for  the  first  mode  of  vibration  in  each 
configuration  are  plotted  as  shown  in  Figure  3.6.  A  listing  of  the  program 
used  to  compute  the  loss  factors  is  in  Appendix  B. 

From  the  data  presented  in  Figure  3.6,  vdscoelastic  thickness  of 
0.38  mm  (0.015  in)  and  constraining  layer  thicknesses  of  3.18  mm  (0.125  in) 
were  selected  for  the  double  layer  configuration.  This  particular 
configuration  estimates  high  damping  over  a  wider  temperature  range 
than  the  other  thickness  combinations. 


TABLE  3.1.  THICKNESSES  USED  IN  CALCULATION  OF  DOUBLE 
LAYER  MODAL  LOSS  FACTORS. 


Configuration  I  2  ^  5 


H2  (mm) 

0.38 

0.38 

0.76 

0.38 

1.14 

HI'  (mm) 

3.18 

2.38 

2.38 

2.38 

2.38 

H2'  (mm) 

0.38 

0.38 

0.76 

1.14 

0.38 

H3'  (mm) 

3.18 

3.18 

2.38 

2.38 

2.38 

The  milled  “floating  element”  plate  uses  the  same  viscoelastic  and 
constraining  layer  thicknesses  as  the  simple  two-layer  configuration.  In  a 
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design  similar  to  that  of  the  pocket  plate,  the  floating  element  and  both 
layers  of  viscoelastic  are  recessed  into  a  milled  opening  as  shown 
previously  in  Figure  3.1  and  further  described  in  the  floating  element 
system  configuration  of  Figure  3.7. 


26 


Single  layer  configuration 


Double  layer  configuration 


Milled  “pocket  plate” 


Plate  with  “floating  element” 


Figure  3.1.  Four  Damping  Treatment  Configurations. 


Figure  3.2.  Elements  of  a  Simple  Sandwich  Damping  System. 
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Temperature  (°F) 


Figure  3.3.  Carpet  Plot  of  Maxirm  m  Loss  Factors  for 
a  Base  Layer  with  HI  =  12.7  mm. 
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CONSTRAINING  LAYER 
VISCOELASTIC 


CONSTRAINING  LAYER 
VISCOELASTIC 


BASE  STRUCTURE 


H3’ 

H2' 


HI' 

H2 


HI 


Figure  3.5.  Configuration  of  the  Double  Constrained 
Layer  System. 
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Estimated  Modal  Loss  Factors 


Temperature  C’F) 


Figure  3.6.  Modal  Loss  Factors  for  Double  Constrained  Layer 
Configurations. 
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IV.  EXPERIMENTAL  RESULTS 


A.  TESTING  ARRANGEMENT 

Experimental  testing  was  performed  on  each  of  the  four  damping 
configurations  and  the  undamped  reference  plate.  In  order  to  approximate 
the  free-free-free-free  boundary  condition,  each  plate  was  suspended  from 
the  roof  of  the  testing  chamber  using  elastic  cords  as  shown  in  Figure  4.1. 
All  of  the  tests  were  performed  in  a  temperature  controlled  environmental 
chamber  which  enabled  temperatures  to  be  maintained  within  ±1  °C.  The 
primary  component  and  user  interface  was  the  Hewlett-Packard  (HP) 
3562A  Dynamic  Signal  Analyzer  (DSA).  The  HP-3562A  was  used  to  provide 
a  swept  sine  signal  to  a  vibration  generator,  and  analyzed  the  returning 
data  signals.  The  HP-3562A  was  used  to  compute  the  frequency  response 
and  coherence  over  a  range  of  50  Hz  to  1050  Hz  using  the  discrete  Fourier 
Transform  in  the  swept  sine  mode.  Ten  averages  were  performed  at  each 
data  point  using  a  frequency  resolution  of  625  mHz  per  step.  The  source 
level  output  to  the  vibration  generator  was  set  at  1.5  volts. 

A  schematic  of  the  testing  apparatus  is  shown  in  Figure  4.2.  Swept 
sine  source  signals  were  fed  from  the  output  jack  of  the  HP-3562A  to  a 
Wilcoxon  F3  vibration  generator  via  the  piezoelectric  output  of  a  Wilcoxon 
PA7C  power  amplifier.  The  vibration  generator  was  mounted  73.48  cm 
(28.93  in)  from  one  end,  and  12.7  cm  (5.0  in)  from  the  front  edge  of  each 
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specimen  as  shown  in  Figure  4.3.  An  integral  force  transducer  was 
mounted  in  the  base  of  the  vibration  generator  to  measure  the  force  input  to 
the  plate.  This  force  signal  was  then  fed  to  input  channel  one  of  the  DSA 
via  a  PCB  462— A  charge  amplifier.  Plate  accelerations  were  recorded  at 
various  points  using  a  PCB  303A-03  accelerometer  as  shown  in  Figure  4.3. 
Acceleration  data  was  fed  to  input  channel  two  of  the  DSA  via  a  PCB  482A05 
power  supply.  Frequency  response  and  coherence  data  was  then  recorded 
on  disk  for  further  analysis. 

Temperatures  within  the  testing  chamber  were  maintained  using  a 
NESLAB  RTE-8  refrigerated  circulating  bath  which  pumped  fluid  through 
a  small  heat  exchanger  in  the  testing  chamber  as  shown  in  Figure  4.2.  In 
order  to  accurately  monitor  the  temperature  of  the  plates,  a  small 
thermocouple  was  inserted  in  the  base  of  each  plate. 


B.  TESTING  PROCEDURE 

1 .  Undamped  reference  plate 

An  undamped,  reference,  frequency  response  measurement 
was  made  at  a  temperature  of  15.6  °C  (60  °F)  to  set  a  standard  response  by 
which  to  measure  the  effectiveness  of  the  damping  treatments.  The 
imdamped  frequency  response  was  recorded  over  a  frequency  range  of 
50  -  1050  Hz  using  a  resolution  of  625  mHz  per  point  in  the  DSA. 
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2.  Damped  plate  measurements 


Frequency  response  measurements  of  the  damped  plates  were 
made  at  a  temperature  of  15.6  °C  (60  °F)  at  several  nodes  on  the  plates  in 
order  to  capture  the  damped  response  of  as  many  modes  as  possible.  A 
representation  of  these  nodes  is  shown  in  Figure  4.3.  Responses  were 
recorded  over  a  frequency  range  of  50  -  1050  Hz  with  a  resolution  of 
625  mHz  per  point  in  the  DSA.  Zoom  measurements  were  also  made  to 
capture  better  data  for  certain  modes.  Modal  loss  factors  were  then 
estimated  from  the  frequency  response  and  coherence  measurements 
using  a  curve-fitting  technique  described  in  Reference  [14]. 


C.  SINGLE  DAMPING  LAYER  RESULTS 

The  single  damping  layer  treatment  was  tested  at  4.44  °C  (40  °F), 

15.6  °C  (60  °F),  and  26.7  °C  (80  °F)  so  that  the  effects  of  temperature  on  the 
damping  treatment  could  be  determined.  A  plot  of  the  single  damping 
layer  frequency  response  at  15.6  °C  is  shown  in  Figure  4.4.  The  single  layer 
damping  treatment  resulted  in  high  damping  with  modal  loss  factors 
ranging  from  0.223  at  53.4  Hz  to  0.091  at  876.6  Hz.  Due  to  the  coupling  of 
modes,  loss  factors  for  all  modes  were  not  measured.  The  frequency 
response  of  the  damped  plate  is  characterized  by  a  frequency  shift  to  the  left 
and  a  smoothing  of  the  frequency  response  when  compared  to  the 
undamped  reference  plate.  The  single  layer  treatment  was  especially 
effective  at  reducing  the  frequency  response  of  a  mode  cluster  between  650 
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and  950  Hz.  The  frequency  band  of  this  cluster  was  shifted  approximately 
200  Hz  with  the  amplitudes  of  the  responses  of  the  modes  being 
dramatically  reduced.  The  single  layer  treatment  was  also  effective  at 
reducing  the  amplitude  of  the  response  peaks  over  the  entire  spectrxim  of 
measurement.  On  average,  the  highest  peaks  of  the  frequency  response  in 
the  undamped  condition  were  reduced  by  25  decibels,  a  reduction  of  17.8 
times. 


The  effect  of  temperature  on  the  damping  was  quite  pronounced  as 
shown  in  Figure  4.5.  As  the  testing  temperature  was  decreased,  the 
viscoelastic  layer  became  stiffer  and  damping  levels  were  increased.  A 
comparative  listing  of  the  loss  factors  at  different  temperature  is  in  Table 
4. 1  and  a  plot  of  the  modal  loss  factors  is  shown  in  Figure  4.6.  Figure  4.5 
shows  the  trend  of  increased  damping  with  temperature  decrease,  and  a 
corresponding  shift  of  modal  frequencies  to  the  right  as  the  viscoelastic 
becomes  stiffer.  These  changes  are  especially  discernible  at  the  lower 
frequencies. 
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TABLE  4.1.  MEASURED  MODAL  LOSS  FACTORS  FOR  THE 


SINGLE  LAYER  AT  DIFFERENT  TEMPERATURE. 


4.44  °C  (40  °F)  15.6  °C  (60  °F)  26.7  °C  (80  °F) 


f(Hz) 

n 

flHzl 

64.3 

0.217 

53.7 

102.8 

0.223 

89.1 

157.4 

0.203 

138.9 

208.4 

0.183 

188.1 

327.5 

0.158 

301.6 

440.1 

0.144 

424.8 

474.9 

0.205 

476.8 

653.1 

0.079 

608.1 

768.0 

0.135 

644.2 

852.0 

0.093 

722.5 

919.0 

0.104 

817.6 

953.5 

0.073 

876.6 

n 

mzi 

n 

0.223 

49.0 

0.117 

0.145 

83.5 

0.089 

0.172 

129.6 

0.082 

0.184 

177.9 

0.072 

0.120 

243.9 

0.062 

0.109 

409.4 

0.053 

0.111 

450.3 

0.066 

0.067 

552.9 

0.068 

0.069 

586.7 

0.052 

0.130 

629.3 

0.042 

0.081 

681.4 

0.096 

0.091 

790.1 

0.049 

858.2 

0.056 

978.1 

0.069 
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D.  DOUBLE  DAMPING  LAYER  RESULTS 


The  double  layer  damping  configuration  was  also  tested  at  4.44  °C, 
15.6  °C,  and  26.7  °C.  The  frequency  response  of  this  configuration  at  15.6  °C 
as  compared  to  the  undamped  reference  plate  is  shown  in  Figure  4.7. 
Damping  in  the  double  layer  configuration  is  also  high,  with  modal  loss 
factors  ranging  form  0.301  at  53.3  Hz  to  0.107  at  832.4  Hz.  Due  to  modal 
coupling  loss  factors  for  all  modes  were  not  measured.  The  frequency 
response  of  the  two-layer  configuration  is  also  characterized  by  a  dramatic 
reduction  in  response  amplitude  and  a  frequency  shift  to  the  left.  The  peak 
undamped  responses  were  reduced  by  an  average  of  27  decibels,  or  a 
reduction  of  22.4  times  from  the  reference  condition. 

The  effect  of  temperature  on  the  double  layer  damping  treatment  is 
shown  in  Figure  4.8.  As  with  the  single  layer  case,  damping  in  the  double 
layer  configuration  increased  with  a  decrease  in  temperature.  This 
configuration  also  shows  the  shift  of  modal  frequencies  to  the  right  as 
temperature  decreases  and  the  viscoelastic  becomes  stiffer.  Modal  loss 
factors  for  the  double  layer  configuration  are  listed  in  Table  4.2,  and  are 
plotted  for  comparison  in  Figure  4.9. 

To  compare  the  effectiveness  of  the  single  layer  and  double  layer 
configurations  their  frequency  responses  at  15.6  °C  are  plotted  in  Figure 
4.10  with  a  plot  comparing  modal  loss  factors  in  Figure  4.11.  The  two 
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responses  are  quite  similar  although  the  double  layer  configuration  does 
show  an  increase  of  approximately  22  percent  in  modal  loss  factor. 


TABLE  4J2.  MEASURED  MODAL  LOSS  FACTORS  FOR  THE 

DOUBLE  LAYER  AT  DIFFERENT  TEMPERATURE 


4.44  =C  (40  °F)  15.6  °C  (60  °F)  26.7  °C  (80  °F) 


flHz) 

n 

mzi 

n 

f(Hz) 

n 

51.9 

0.273 

55.3 

0.301 

49.7 

0.202 

122.5 

0.224 

87.1 

0.215 

81.9 

0.188 

212.7 

0.187 

142.1 

0.217 

127.6 

0.156 

290.5 

0.197 

190.5 

0.212 

172.4 

0.117 

382.7 

0.174 

297.8 

0.139 

278.2 

0.097 

494.2 

0.144 

366.8 

0.154 

317.7 

0.077 

558.2 

0.169 

419.3 

0.125 

3%.l 

0.070 

602.2 

0.198 

441.2 

0.100 

428.5 

0.077 

642.4 

0.157 

618.7 

0.098 

602.0 

0.067 

752.8 

0.168 

680.9 

0.096 

638.4 

0.072 

815.4 

0.159 

715.5 

0.060 

738.0 

0.071 

896.3 

0.156 

832.4 

0.107 

802.6 

0.077 

968.3 

0.130 

846.1 

0.050 
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E.  POCKET  PLATE  RESLXTS 


The  milled  pocket  plate  was  constructed  as  previously  shown  in 
Figure  3.3.  The  constraining  layer,  or  cover  plate,  was  welded  in  place 
using  tack  welds  in  an  attempt  to  keep  the  damping  material  away  from  the 
heat  of  welding,  and  the  cover  plate  from  warping,  instead  of  using  a 
continuous  weld  bead  as  was  previously  attempted  [Ref.  2].  The  cover  plate 
was  welded  to  the  base  at  the  corners,  at  the  midpoint  of  the  short  side  and 
at  three  equally  spaced  locations  along  the  long  dimension  as  shown  in 
Figure  4.12.  Following  welding  the  plate  was  tested  to  ensure  that  the 
viscoelastic  had  not  been  damaged  by  the  heat  of  welding. 

The  pocket  plate  was  tested  at  15.6  °C  (60  °F)  and  the  frequency 
response  is  shown  in  Figure  4.13.  The  response  indicates  that  the 
viscoelastic  laye"  was  not  damaged  by  welding  and  that  good  damping  was 
attained.  Modal  loss  factors  ranged  from  0.067  at  62.1  Hz  to  0.090  at  923  Hz. 
Although  damping  is  good,  it  is  approximately  half  that  of  the  single  layer 
configuration.  One  reason  for  this  is  that  the  viscoelastic  material  does  not 
completely  cover  the  base  structure.  Another  reason  is  the  presence  of  the 
welded  cover  plate.  Due  to  the  welded  conditions  the  cover  plate  cannot 
induce  shear  deformation  in  the  viscoelastic  layer  as  well  as  a  true 
constraining  layer,  and  therefore  produces  less  damping  than  the  single 
layer  configuration.  The  effects  of  the  welded  cover  plate  are  especially  felt 
in  modes  below  300  Hz  where  the  frequency  response  is  quite  peaked.  The 
response  curve  becomes  more  rounded  and  the  effects  of  the  damping  layer 
are  seen  as  frequency  increases. 
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The  modal  loss  factors  for  the  pocket  plate  are  listed  in  Table  4.3  and 
are  plotted  in  Figure  4.14.  Even  though  the  damping  is  less  than  the  single 
layer,  the  plate  is  still  adequately  damped  as  shown  in  the  frequency 
response  plot  in  Figure  4.13.  In  this  configuration  the  modal  loss  factors 
remained  relatively  constant  throughout  the  testing  spectrum.  The 
increase  in  modal  loss  factor  values  above  800  Hz  is  due  primarily  to  modal 
coupling,  and  the  measured  modal  loss  factors  in  this  range  are  not 
reliable. 


TABLE  4.3  IVIEASURED  MODAL  LOSS  FACTORS  FOR  THE 
POCKET  PLATE  AT  15.6°C. 


f(Hz) 

n 

f(Hz) 

n 

62.1 

0.067 

565.3 

0.043 

93.4 

0.042 

625.7 

0.047 

142.0 

0.056 

643.2 

0.049 

194.9 

0.060 

687.0 

0.044 

308.9 

0.051 

841.9 

0.056 

437.5 

0.041 

891.5 

0.077 

495.6 

0.036 

923.0 

0.090 
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F. 


FLOATING  ELEMENT  RESLT.TS 


The  milled  “floating  element”  configuration  was  constructed  as 
previously  shown  in  Figure  3.7.  The  cover  plate  was  welded  in  a  fashion 
similar  to  the  pocket  plate  as  shown  in  Figure  4.12.  The  center 
constraining  layer,  or  floating  element,  was  made  slightly  smaller  than  the 
surrounding  structure  thus  allowing  the  floating  element  to  act  as  a  “true” 
constraining  layer. 

The  frequency  response  of  the  floating  element  configuration  at 
15.6  °C  (60  °F)  is  shown  in  Figure  4.15.  The  floating  element  is  quite 
effective  as  the  response  shows  a  good  reduction  in  peak  modal  response. 
Measured  modal  loss  factors  range  from  0  089  at  66  Hz  to  0.064  at  935  Hz.  A 
listing  of  measured  modal  loss  factors  is  in  Table  4.4  and  are  plotted  in 
Figure  4.16.  As  with  the  previous  cases,  the  frequency  response  of  the 
floating  element  configuration  is  characterized  by  a  frequency  shift  to  the 
left  and  a  smoothing  of  the  response  as  frequency  increases. 

In  a  comparison  of  the  pocket  plate  and  floating  element 
configurations,  the  frequency  responses  are  plotted  in  Figure  4.17.  A 
comparison  of  modal  loss  factors  for  the  two  configurations  is  shown  in 
Figure  4.18.  The  two  frequency  response  plots  are  similar,  however,  the 
frequency  response  of  the  floating  element  configuration  is  more  rounded 
than  that  of  the  pocket  plate.  The  major  difference  between  the  two 
configurations  is  seen  in  Figure  4.18.  Modal  loss  factors  for  the  floating 
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element  show  an  average  increase  of  25  percent  over  those  of  the  pocket 
plate.  Reasons  for  this  increase  are  the  added  constraining  effect  of  the 
floating  element  and  additional  layer  of  damping  material  present. 


TABLE  4.4  MEASURED  MODAL  LOSS  FACTORS  FOR  THE 
FLOATING  ELEMENT  AT  15.6°C. 


f(Hz) 

n 

f(Hz) 

n 

66.0 

0.089 

508.1 

0.089 

104.9 

0.040 

540.6 

0.115 

154.0 

0.094 

618.7 

0.121 

211.4 

0.058 

656.7 

0.100 

272.8 

0.063 

708.6 

0.088 

323.5 

0.075 

857.6 

0.089 

435.8 

0.090 

935.4 

0.064 
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Fi^re  4.1.  Testing  Configuration  in  Testing  Chamber. 
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Figure  4.2.  Schematic  Diagram  of  Testing  System. 
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REPRESENTATIVE  ’ 
MEASUREMENT  POINTS 


Figure  4.3.  Shaker  and  Accelerometer  Locations. 
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Figure  4.5.  Frequency  Response  of  the  Single  Layer 

Configuration  at  Different  Temperatures. 

[ - :  4.44  °C, - :  15.6  °C, - :  26.7  °C] 


8 


Modal  Ix)ss  Factors 


Frequency  (Hz) 


Figure  4.6.  Comparison  of  Modal  Loss  Factors  for  the  Single 
Layer  Configuration  at  Different  Temperatures. 
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\ _ ! _ 1 _ 1 _ 1 _ ^ _ I _ I _ \ _ \ _ 

FxdXY  50  HZ  DOUBLE  LAYER  15.6  C  1.05k 


Figure  4.7.  Frequency  Response  of  the  Double  Layer  at  15.6  °C. 
[ - :  undamped  response  , - :  damped  response] 
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I _ I _ I _ I _ L_ _ I _ 1 _ I _ i _ I  I 

FxdXY  5C  Hz  4.44  /  15.5  /  26.7  C  1.05k 


Figure  4.8.  Frequency  Response  Comparison  of  the  Double 
Layer  Configuration  at  Different  Temperatures. 

[ - :  4.44  °C, - :  15.6  °C  , - :  26.7  °C] 
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Modal  Loss  Factors 
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Figure  4.9.  Comparison  of  Modal  Loss  Factors  for  the  Double 
Layer  at  Different  Temperatures. 
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FxdXY  50  Hz  SINGLE  /  DOUBLE  1.05k 


Figure  4.10.  Frequency  Response  of  the  Double  and 
Single  Layer  at  15.6  °C. 
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Frequency  (Hz) 

Figure  4.11.  Modal  Loss  Factors  for  Single  and  double  layer 
Configurations. 


Figure  4.12.  Location  of  Tack  Welds  on  the  Cover  Plate 
of  the  Pocket  Plate  Configuration. 
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Fxd  Y  50  HZ  POCKET  PLATE  15.6  C  1.05k 


Figure  4.13.  Frequency  Response  of  the  Pocket  Plate 
Configuration  at  15.6  °C. 

[ - :  reference  plate  , - :  pocket  plate  ] 
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Fxd  Y  50  HZ  FLOATING  ELMT  15.6  C  1 . OBk 


Figure  4.15.  Frequen^'y  Response  of  the  Floating  Element 
at  15.6  °C  . 

[ - ;  undamped  plate  , - :  floating  element  ] 
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Modal  Ivoss  Factors 
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Figure  4.16.  Modal  Loss  Factors  for  the  Floating  Element 
Configuration. 
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Fxd  Y  50  Hz  POCKET  /  FLOATING  1.05k 


Figure  4.17.  Comparison  of  Frequency  Responses  for  the 

Pocket  Plate  and  Floating  Element  Configurations. 
[ - :  pocket  plate  , - :  floating  element  ] 
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Modal  Ix)s.s  Factors 
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Figure  4.18.  Modal  Loss  Factors  for  the  Pocket  Plate 
and  Floating  Element  Configurations. 
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V.  FINITE  ELEMENT  RESULTS 


A .  UNDAMPED  REFERENCE  PLATE 

The  first  step  in  the  finite  element  analysis  procedure  was  to  model 
and  analyze  the  undamped  reference  plate  for  its  modal  frequencies  and 
frequency  response.  The  finite  element  model  was  generated  using 
PATRAN,  a  computer  aided  interactive  graphics  program  developed  by 
PDA  Engineering.  PATRAN  is  widely  used  for  conceptual,  preliminary, 
and  detailed  design  and  analysis  of  complex  systems.  One  of  PATRAN’s 
major  advantages  is  in  the  interactive  construction  of  finite  element  models 
for  use  by  MSC/NASTRAN,  and  its  ability  to  display  MSC/NASTRAN 
results  in  an  easily  understood  graphic  format  [Ref.  15]. 

The  reference  plate  was  modeled  using  84  plate  (QUAD4)  elements  as 
shown  in  Figure  5.1.  The  QUAD4  element  is  an  isoparametric  element 
with  four  nodes,  one  at  each  corner  of  the  element  [Ref.  7].  A  normal  mode 
extraction  was  then  performed  in  order  to  compare  numerical  results  with 
experimentally  determined  modal  frequencies.  Once  satisfactory 
agreement  between  the  modal  frequencies  calculated  in  NASTRAN  and 
those  obtained  experimentally  was  attained,  a  direct  frequency  response 
calculation  was  performed  in  NASTRAN.  This  frequency  response  was 
then  used  as  the  reference  response  for  further  finite  element  models 
incorporating  viscoelastic  damping  treatments.  The  excitation  for  the 


frequency  response  calculation  was  a  sinusoidal  force  with  an  amplitude  of 
1.0  applied  at  node  66.  The  response  point  was  node  58  as  shown  in  Figure 
5.1.  These  two  nodes  correspond  to  the  points  on  the  plates  used  in  the 
experimental  portion  of  the  research  where  the  vibration  generator  and 
accelerometer  were  attached. 

R  SINGLE  DAMPING  LAYER 

The  modeling  of  the  single  constrained  layer  damping  system  was 
done  using  techniques  described  by  Johnson  and  Kienholz  [Ref  1].  As 
shown  in  Figure  5.2,  the  viscoelastic  layer  is  modeled  using  solid  (HEXA) 
elements,  while  the  base  layer  and  constraining  layer  were  modeled  using 
QUAD4  elements.  The  HEXA  element  is  a  solid,  isoparametric  element 
having  eight  nodes,  one  at  each  corner  of  the  element  with  three 
translational  degrees  of  freedom  at  each  node  [Ref.  7].  The  use  of  solid 
elements  for  the  viscoelastic  layer  allows  the  strain  energy  due  to  shearing 
to  be  adequately  represented.  Plate  elements  are  used  in  the  base  layer  and 
constraining  layer  because  of  their  ability  to  account  for  stretching  and 
bending  deformations.  The  plate  element  allows  its  nodes  to  be  offset  from 
the  plate’s  center  to  the  surface  of  the  plate,  coincident  with  the  corner 
nodes  of  the  solid  viscoelastic  elements  [Ref.  1].  Thus,  the  single 
constrained  layer  system  was  modeled  using  only  two  layers  of  nodes,  a 
simple  process  in  PATRAN.  For  the  single  layer  configuration,  a  model 
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having  84  elements  per  layer,  and  an  element  meshing  scheme  shown  in 
Figure  5.2  was  used. 

Once  the  damped  plate  had  been  modeled,  normal  mode  extractions 
were  made  using  MSC/NASTRAN.  Five  separate  runs  were  conducted 
using  the  material  properties  of  ISD-112  at  50,  200,  500,  800,  and  1000  Hz 
and  at  a  temperature  of  15.6  °C  (60  °F).  In  addition  to  the  modal 
frequencies,  the  strain  energy  in  the  viscoelastic  elements  and  the  entire 
model  were  output  from  NASTRAN. 

Since  the  shear  modulus  of  a  viscoelastic  material  changes  with 
frequency,  it  was  necessary  to  estimate  the  actual  modal  frequencies  of  the 
damped  plate  using  an  interpolation  procedure  outlined  by  Johnson  and 
Kienholz  [Ref.  16].  The  first  step  of  the  interpolation  process  was  to  plot  the 
shear  modulus  of  ISD-112  versus  frequency  from  5  to  1000  Hz,  Then,  for 
the  first  mode,  using  NASTRAN  results  based  on  ISD-112  material 
properties  a  50  Hz,  the  first  modal  frequency  predicted  by  NASTRAN  and 
the  corresponding  shear  modulus  were  plotted.  The  same  was  then  done 
using  the  first  natural  frequency  predicted  by  normal  mode  extraction 
based  on  ISD-112  material  properties  at  frequencies  of  500  and  1000  Hz.  A 
curve  was  then  passed  through  these  three  points.  The  point  where  the 
NASTRAN  modal  frequencies  for  the  first  mode  intersected  the  ISD-112 
shear  modulus  curve  was  taken  to  be  the  interpolated  modal  frequency  of 
the  single  damping  layer  configuration.  A  plot  of  the  intersection  of  these 
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two  curves  for  the  first  and  second  modes  is  shown  in  Figure  5.3.  This 
interpolation  process  was  then  repeated  for  each  mode  through  1000  Hz. 

Once  the  interpolated  modal  frequencies  were  found,  the  modal 
strain  energy  equations  (2.12)  and  (2.13)  were  used  to  compute  modal  loss 
factors  for  the  single  layer  configuration.  A  set  of  modal  loss  factors  was 
computed  based  on  the  modal  strain  energies  computed  using  viscoelastic 
properties  at  reference  frequencies  of  50  ,200,  500,  and  800  Hz.  A  set  of 
composite  modal  loss  factors  for  the  modal  frequencies  near  these  reference 
frequencies  was  then  selected.  The  resulting  modal  loss  factors  are  shown 
in  Table  5.1  and  are  plotted  versus  frequency  in  Figure  5.4.  As  seen  in 
Figure  5.4,  the  modal  strain  energy  method  is  predicting  high  damping  for 
this  configuration  with  an  average  modal  loss  factor  of  0.195. 


TABLE  5.1.  ESTIMATED  MODAL  LOSS  FACTORS  FOR  THE 
SINGLE  LAYER  USE^G  THE  MODAL  STRAIN 
ENERGY  METHOD 


f^Hz) 

n 

f(Hz) 

a 

63 

0.249 

517 

0.237 

95 

0.181 

586 

0.217 

153 

0.234 

632 

0.164 

202 

0.212 

663 

0.155 

285 

0.229 

712 

0.180 

322 

0.214 

827 

0.134 

447 

0.200 

869 

0.152 

463 

0.108 

883 

0.115 

483 

0.258 

65 


Using  the  set  of  composite  modal  loss  factors,  the  modal  frequency 
response  of  the  damped  plate  was  computed  using  MSC/NASTRAN.  Modal 
damping  was  introduced  to  the  model  using  the  SDAMP  option  in  the  Case 
Control  Deck  and  the  TABDMPl  damping  table  in  the  Bulk  Data  Deck  as 
described  in  the  MSC/NASTRAN  Handbook  for  Dynamic  Analysis  [Ref.  17]. 
Since  NASTRAN  uses  a  linear  interpolation  between  points  in  the  damping 
table  to  describe  the  modal  damping  in  the  model  [Ref  7],  a  simple  curve  fit 
was  applied  to  the  set  of  composite  modal  loss  factors  as  shown  in 
Figure  5.4.  Points  from  this  curve  fit  were  then  used  in  the  NASTRAN 
damping  table.  To  compute  the  modal  frequency  response,  a  unit  excitation 
force  was  applied  at  the  same  node  as  the  undamped  plate,  and  the  node 
used  for  the  response  was  also  the  same  as  the  undamped  plate. 

The  results  of  the  modal  frequency  response  calculations  are  shown 
in  Figure  5.5.  The  dashed  line  represents  the  undamped  reference  plate, 
and  the  solid  lii.e  represents  the  modal  frequency  response  of  the  single 
layer  configuration.  Material  properties  at  200  Hz  were  used  for  the 
ISD-112  damping  material.  The  first  thirty  modes  were  used  in  the  modal 
summation  for  the  response.  A  listing  of  the  MSC/NASTRAN  data  deck 
used  to  compute  the  modal  frequency  response  is  in  Appendix  D. 

The  modal  loss  factors  estimated  using  the  modal  strain  energy 
method  are  compared  to  those  measured  experimentally  for  the  single  layer 
configuration  in  Figure  5.6.  The  estimated  loss  factors  are  greater  than  the 


66 


experimentally  determined  loss  factors  throughout  the  spectrum  of 
interest,  and  especially  in  the  middle  frequencies. 

The  modal  frequency  response  of  the  single  layer  configuration  is 
compared  to  the  experimentally  measured  frequency  response  in 
Figure  5.7.  The  comparison  was  accomplished  by  normalizing  both  the 
experimentally  determined  frequency  response  and  the  frequency  response 

computed  in  NASTRAN.  Both  responses  were  normalized  using  a  value  of 
in/s0C^ 

1.0  — .  The  effects  of  the  greater  loss  factors  estimated  by  the  modal 

strain  energy  method  are  obvious  as  the  level  of  the  predicted  response  is 
lower  than  the  measured  response.  The  shift  in  frequency  between  the  two 
curv'es  is  due  to  the  finite  element  model  being  inherently  stiffer  than  the 
actual  system.  The  correlation  between  the  two  curves  is  especially  good 
below  250  Hz  as  this  is  where  the  differences  between  estimated  and 
measured  modal  loss  factors  are  the  least. 


C .  DOUBLE  DAMPING  LAYER 


The  modeling  of  the  double  constrained  layer  damping  system  was 
accomplished  as  shown  in  Figure  5.8.  The  double  layer  configuration 
consists  of  a  base  layer  modeled  with  offset  QUAD4  elements,  two 
viscoelastic  layers  consisting  of  HEXA  elements,  and  the  top  constraining 
layer  modeled  with  offset  QUAD4  elements.  The  middle  constraining  layer 
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was  modeled  using  three  layers  of  HEXA  elements  in  order  to  give  this 
layer  the  stiffness  necessary  to  act  as  a  constraining  layer.  The  model  was 
meshed  using  60  elements  in  each  layer  as  shown  in  Figure  5.8. 

Using  this  model,  the  modal  strain  energy  method  was  applied  to 
determine  approximate  modal  frequencies  and  loss  factors  for  the  double 
layer  configuration.  To  determine  the  loss  factors,  normal  mode 
extractions  were  performed  using  reference  frequencies  of  50,  200,  500,  800, 
and  1000  Hz.  A  composite  set  of  modal  loss  factors  for  the  double  layer 
system  was  then  compiled  based  on  the  estimated  modal  frequency’s 
relation  to  the  reference  frequency  used  to  calculate  modal  strain  energies. 
This  composite  set  of  modal  loss  factors  is  listed  in  Table  5.2,  and  is  plotted 
in  Figure  5.19  as  a  comparison  to  the  experimentally  determined  modal 
loss  factors  for  the  double  layer  configuration.  The  estimated  modal  loss 
factors  for  the  double  layer  show  high  damping,  but  they  compare  favorably 
with  those  measured  experimentally. 

The  modal  frequency  response  of  the  double  layer  configuration  was 
computed  in  a  manner  similar  to  the  single  layer  in  that  smoothed  loss 
factor  data  was  used  in  the  MSC/NASTRAN  damping  table.  Likewise,  a 
unit  excitation  force  was  applied,  and  the  first  thirty  modes  were  used  in 
the  modal  summation.  The  results  of  the  modal  frequency  response 
calculation  are  shown  in  Figure  5.10.  The  dashed  line  represents  the 
undamped  response,  and  the  solid  line  represents  the  frequency  response  of 
the  double  layer  configuration. 
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TABLE  5.2.  ESTIMATED  MODAL  LOSS  FACTORS  FOR  THE 


DOUBLE  LAYER  USING  THE  MODAL  STRAIN 
ENERGY  MEIHOD. 


ffHz) 

n 

f^Hz) 

n 

59 

0.278 

544 

0.209 

89.5 

0.219 

578 

0.141 

140 

0.239 

593 

0.131 

189 

0.213 

659 

0.177 

258 

0.214 

754 

0.102 

300 

0.205 

783 

0.076 

402 

0.179 

802 

0.144 

427 

0.176 

956 

0.077 

456 

0.264 

969 

0.111 

486 

0.242 

994 

0.048 

The  frequency  response  calculated  using  NASTRAN  was  compared 
to  the  experimentally  determined  frequency  response  of  the  double  layer 
configuration  as  shown  in  Figure  5.11.  Comparison  of  the  two  frequency 
response  curves  shows  similarity  in  form,  but  a  much  lower  response  level 
for  the  numerically  determined  response.  Once  again  this  could  be  due  to 
the  higher  damping  predicted  by  the  modal  strain  energy  method  and  the 
inherently  higher  stiffness  of  the  finite  element  model. 
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D.  POCKET  PLATE  RESLXTS 


The  pocket  plate  configuration  was  modeled  with  the  same  offset 
plate  elements  and  solid  viscoelastic  elements  as  the  single  layer 
configuration.  However,  the  pocket  plate  required  the  modeling  of  the 
milled  structure  around  the  viscoelastic  material  and  the  welds  between 
the  cover  plate  and  milled  plate.  A  representation  of  the  model  is  shown  in 
Figure  5.12.  The  base  structure,  cover  plate,  and  the  structure  immediately 
around  the  cover  plate  was  modeled  using  offset  plate  elements.  The 
vdscoclastic  material  and  the  portion  of  the  structure  immediately  adjacent 
to  it  were  modeled  using  the  solid  HEXA  elements.  Since  the  viscoelastic 
material  and  cover  plate  are  physically  separated  from  the  surrounding 
plate,  except  where  the  viscoelastic  is  adhered  to  the  base  structure,  care 
was  necessary  in  creating  the  finite  element  mesh. 

The  model  was  created  in  PATRAN  using  PATRAN’s  node  editing 
and  equivalencing  capabilities  [Ref.  18].  This  allowed  the  generation  of  a 
finite  element  mesh  with  two  nodes  at  the  same  geometric  point  in  space. 
Using  this  node  editing  technique,  a  mesh  was  created  which  allowed  the 
viscoelastic  and  cover  plate  to  vibrate  separately  from  the  surrounding 
structure,  yet  at  the  same  time,  keep  the  number  of  elements  and  nodes  in 
the  model  to  a  minimum.  The  welded  points  on  the  cover  plate  were  also 
modeled  using  node  editing  techniques.  At  weld  points  the  finite  element 
node  on  the  cover  plate  was  equivalenced  with  its  corresponding  node  on  the 
base  structure,  resulting  in  a  single  node  and  a  connection  between  an 
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otherwise  separate  base  structure  and  cover  plate.  At  non— welded  points 
on  the  cover  plate  there  were  two  nodes  at  the  same  geometric  point;  one  to 
represent  the  cover  plate,  and  the  other  to  represent  the  base  structure.  The 
model  was  meshed  using  a  5x11  mesh  resulting  in  40  elements  in  the  cover 
plate,  viscoelastic  and  base  layer  as  shown  in  Figure  5.12. 

Using  this  modeling  scheme,  the  modal  strain  energy  method  was 
employed  to  estimate  the  modal  frequencies  and  loss  factors.  Normal  mode 
extractions  were  made  using  viscoelastic  material  properties  at  50,  200,  500, 
800,  and  1000  Hz.  Using  these  reference  frequencies  a  composite  set  of 
modal  loss  factors  was  obtained.  These  loss  factors  are  listed  in  Table  5.3 
and  are  plotted  versus  frequency  in  Figure  5.13.  The  estimated  loss  factors 
give  good  damping  over  the  spectrum  of  interest  with  an  average  modal  loss 
factor  of  0.075.  As  with  the  pre\dous  cases,  damping  values  used  for  the 
modal  frequency  response  calculation  came  from  a  curve-fit  to  the  set  of 
composite  modal  loss  factors. 

The  modal  frequency  response  of  the  pocket  plate  was  computed 
using  the  first  30  modes  and  viscoelastic  material  properties  at  200  Hz.  The 
resulting  estimated  frequency  response  is  shown  in  Figure  5.14.  The 
response  shows  a  definite  frequency  shift  to  the  left  along  with  good 
damping  of  the  frequency  response  when  compared  to  the  undamped 
response. 
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The  estimated  modal  loss  factors  and  modal  frequency  response  for 
the  pocket  plate  were  compared  to  those  measured  experimentally.  The 
loss  factor  comparison  is  shown  is  Figure  5.15  and  the  frequency  response 
comparison  is  shown  in  Figure  5.16.  The  estimated  modal  loss  factors  are 
higher  than  those  measured  experimentally,  however,  the  frequency 
responses  compare  quite  favorably  with  each  other.  The  frequency 
response  curv'e  calculated  through  finite  element  analysis  has  a  lower 
response  level  and  a  frequency  shift  to  the  right  of  the  measured  frequency 
response.  This  is  expected  due  to  the  increase  in  damping  predicted  by  the 
modal  strain  energy  method  and  by  the  fact  that  the  finite  element  model  is 
inherently  stiffer  than  the  physical  system 


TABLE  5.3.  ESTIMATED  MODAL  LOSS  FACTORS  FOR  THE 
POCKET  PLATE  USING  THE  MODAL  STRAIN 
ENERGY  METHOD. 


f(Hz') 

n 

f(Hz) 

n 

67 

0.046 

511 

0.113 

109 

0.034 

588 

0.104 

162 

0.091 

625 

0.086 

217.5 

0.063 

668 

0.089 

295 

0.091 

713 

0.077 

333 

0.047 

831 

0.077 

453 

0.07 

834 

0.065 

477 

0.093 

868 

0.051 
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Figure  5.1.  Finite  Element  Model  of  the  Undamped 
Reference  Plate. 
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EXCITATION  POINT 
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Figure  5.3.  Interpolation  of  the  First  and  Second  Modal  Frequencies 


for  the  Single  Constrained  Layer  Configuration. 
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Figure  5.4.  Estimated  Modal  Loss  Factors  for  the  Single  Layer 
Configuration  with  the  Curve  Fit  Used  for  the 
MSC/NASTRAN  Damping  Table. 
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Figure  5.5.  Calculated  Modal  Frequency  Response  of  the  Single 
Layer  Configuration  Using  NASTRAN. 

[ - :  reference  plate  , - ;  single  layer] 
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Figure  5.6.  Comparison  of  Estimated  and  Measured  Modal  Loss 
Factors  for  the  Single  Layer  Configuration. 


78 


FREQUENCY  RESPONSE  (DB) 


Figure  5.7.  Comparison  of  Estimated  and  Measured  Frequency 
Response  for  the  Single  Layer  Configuration. 
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Figure  5.9.  Modal  Loss  Factors  for  the  Double  Layer  Configuration 
as  Determined  from  the  Modal  Strain  Energy  Method 
and  Determined  Experimentally. 
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ACCELERATIOf^l  (Db) 


Figure  5.10.  Calculated  Modal  Frequency  Response  for  the  Double 
Layer  Configuration  Using  NASTRAN. 

[ - :  reference  plate  , - :  double  layer] 
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Figure  5.11.  Comparison  of  the  Experimentally  Determined 
and  Numerically  Predicted  Frequency  Rcr.ponses 
for  the  Double  Layer  Configuration. 


83 


.  :  r  - ; — - 1 . . ; - 'T' — ^ - r 

100  200  500  400  500  600  700  800  900  1000 

FREQUENCY  (HZ) 


Figure  5.14.  Estimated  Frequency  Response 

for  the  Pocket  Plate  Configuration. 

[ - :  reference  plate  , - :  pocket  plate] 
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Figui'e  5.15  Comparison  of  Experimentally  Measured  and 
Estimated  Modal  Loss  Factors  for  the  Pocket 
Plate  Configuration. 
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FREQUENCY  RESPONSE  (DB) 


Figure  5.16.  Comparison  of  Experimental  and  Predicted 
Frequency  Responses  for  the  Pocket  Plate 
Configuration. 
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VL  CONCLUSIONS 


Constrained  viscoelastic  layer  damping  is  an  extremely  effective 
method  for  reducing  broadband  vibration.  In  each  of  the  experimental 
cases  the  peak  amplitudes  of  frequency  response  were  reduced  by 
approximately  25  decibels,  a  reduction  of  18  times  below  the  undamped 
reference  plate. 

In  a  comparison  of  experimentally  determined  modal  loss  factors  for 
the  four  treatments,  the  double  layer  configuration  yields  the  largest 
damping  and  the  pocket  plate  the  least  as  shown  in  Figure  6.1.  Of 
particular  note  are  the  performance  of  the  pocket  plate  and  floating  element 
configurations.  Although  they  are  not  ideal  configurations  in  terms  of 
“true”  constrained  viscoelastic  layer  damping,  the  damping  levels  achieved 
are  quite  satisfactory.  As  was  previously  reported  in  Section  IV,  the 
floating  element  configuration  yielded  an  average  increase  of  25  percent 
over  the  modal  loss  factors  of  the  pocket  plate.  It  is  also  noted  that  the 
average  modal  loss  factor  for  the  pocket  plate  is  approximately  50  percent  of 
the  average  modal  loss  factor  of  the  single  layer  configuration.  Similarly, 
the  average  modal  loss  factor  for  the  floating  element  configuration  is 
approximately  50  percent  that  of  the  double  layer  treatment. 

The  modal  strain  energy  method  tends  to  overpredict  the  modal  loss 
factors  by  as  much  as  50  percent  for  these  highly  damped,  thick  plates. 
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Figure  6.2  shows  that,  although  the  modal  strain  energy  method  predicted 
modal  loss  factors  greater  than  those  measured,  the  same  relative 
differences  in  modal  loss  factor  between  damping  configurations  are 
maintained.  The  estimated  modal  loss  factors  for  the  pocket  plate  are 
approximately  60  percent  less  than  those  of  the  single  layer  configuration. 
This  indicates  that  there  is  a  consistency  between  the  damping  values 
predicted  by  the  modal  strain  energy  method  and  those  of  the  actual 
physical  system. 

There  are  several  possible  reasons  for  the  differences  between  the 
experimentally  determined  modal  loss  factors  and  those  estimated  by  the 
modal  strain  energ>’  method.  The  first  is  that  the  material  properties  of 
ISD-112  as  reported  by  the  3M  Corporation  on  the  reduced  frequency 
nomogram  may  not  be  consistent  with  the  material  properties  actually 
present  in  the  material  used.  Since  the  numerical  analysis  was  based  on 
the  reported  material  properties  this  is  a  possible  source  of  uncertainty  in 
the  results. 

The  second  source  of  uncertainty  is  in  the  adhesion  of  the  ISD-112  to 
the  aluminum  plates.  Although  the  plates  were  clean  when  the  ISD-112 
was  applied,  it  was  noted  that  the  adhesive  qualities  of  the  ISD— 112  were  not 
uniform  throughout  the  material.  A  lack  of  adhesion  may  cause  a 
reduction  in  the  damping  capability  of  the  system. 
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Another  source  for  the  difference  between  the  experimental  results 
and  numerical  results  lies  within  the  finite  element  model.  The  number  of 
elements  used  in  the  model  greatly  affects  its  “stiffness.”  As  the  number  of 
elements  is  increased  the  model  should  become  less  stiff  and  results  are 
expected  to  approach  those  that  are  measured.  Also,  the  type  of  element 
used  to  model  the  base  layer  and  constraining  layer  may  have  an  effect.  In 
this  research  plate  elements  were  used  to  model  both  the  base  and 
constraining  layers.  It  is  possible  that  one  or  more  layers  of  solid  elements 
may  produce  results  that  agree  better  with  experimental  results. 

One  drawback  to  the  modal  strain  energy  method  in  design  is  the 
large  amount  of  CPU  time  required  for  normal  mode  extraction  and  modal 
frequency  response;  especially  in  complex  structures  with  a  large  number 
of  elements.  Therefore,  altnough  the  modal  strain  energy  method  is  good 
for  analyzing  a  design,  it  may  face  a  big  difficulty  to  be  used  for  design 
optimization  due  to  the  large  amount  of  CPU  time  required  for  normal 
mode  extraction  and  modal  frequency  response  calculations. 


31 


Modal  I  vOss  Factor 


Figure  6. 
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Figure  6.2.  Comparison  of  Estimated  Modal  Loss  Factors  for 
the  Single  Layer,  Double  Layer,  and  Pocket  Plate 
Configurations. 


93 


\1L  RECOMIVIENDATIONS 


The  thick  plate  used  in  this  research  is  a  generic  model  of  many 
physical  systems  that  may  see  use  in  naval  application.  There  are  several 
areas  which  deserve  more  research  and  clarification,  including  the 
following: 

•  Model  the  floating  element  configuration  in  finite  elements  and 
check  the  effectiveness  of  the  modal  strain  energy  method  in 
predicting  modal  loss  factors  for  this  configuration. 

•  Investigate  the  relationship  between  mode  shapes  and  damping 
values.  It  seems  some  modal  damping  values  are  very  low  due  to 
their  twisLlng  mode  shapes. 

•  For  pocket  plate  and  floating  element  configurations,  weld  the  w’hole 
cover  plate  (continuous  weld)  and  investigate  the  damping 
characteristics  to  compare  the  vibration  reduction  capabilities  of 
the  two  different  milled  plate  damping  treatments. 

•  Investigate  methods  for  improving  the  adhesion  of  ISD-112  to  the 
base  structure  and  constraining  layer. 
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APPENDIX  A 


FORTRAN  PROGRAM  USED  TO  COMPILE  MODAL  LOSS  FACTORS 
FOR  THE  SINGLE  DAMPING  LAYER  DESIGN 

This  program  used  the  Ross-Kerwin-Ungar  equations  of  Section  III 
to  compute  estimated  modal  frequencies  and  loss  factors  for  the  single  layer 
configuration.  Modal  frequencies  for  an  undamped  plate  are  read  from  a 
data  file  and  estimated  modal  frequencies  and  loss  factors  for  various  layer 
thicknesses  are  output  to  another  file.  Material  properties  of  ISD-112  are 
computed  using  University  of  Dayton  data  and  curve-fitting  equations  to  the 
reduced  frequency  nomogram  [Ref.  9,12].  The  units  used  in  this  program 
are  pounds,  inches,  and  seconds. 
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rrair.i'.Aii  (Atari 

c  mil  croiiia'.fi  la  m  cAic  iuAir  r.fGirn  i  oa-r  (  ai  ihh:;  i  dh  vAniour, 

c  1 1 1 II' u; A I  cia  a  amd  iuuh  a  iii  a  CdiiT)  i  ka  j  hi  d  i  a"!  i  li  v  i  ol  iii  i  Aa  1 1  c 

(  HMii'iMi,  5t:)iiM  ir.iiK.  iiii.  Kiia''.  I. riciiiii  uiK'-Aii  LiJUAiiDiia.  am.irn  idiS 
I'  i.'.iiHi:',  Ai:i  (lir.niiiu  i  hk  a  i  i  mi'Lka  nii;t  iiAiic.r.  ui'  :u  loo  ntoniis 
I  lAiaiiiiiii.  iiAar  i'lak  iii!CKm.”.'ji: 3  vaiiv  mcimi  n.;.?*;  10  n  75  iiicnra 
c  111  u  1  iiCKi  niai  I  3 .  MiK  rACH  nA3£  riAir  111KK111.33.  nit  vir.aoi  la3i  ic 
c  I II 1  fiaiL  33  is  vAu'icn  ri-’on  u.ois  10  o.iioo  11111113  lii  ii  0 1  3  iiaii 

L  I  luaiLiinii  3  .  f  UK  LAcii  VI  scdci  AS  1  ic  iHitKin53  me  cons  I  RA I  III  Ik; 

c  LAUK  micaair;.:,  is  varilo  riioii  o.oajs  10  u.cs  iiitiiLS  111 
c  0  0623  men  iiicRciirins. 

C  mis  rilOoRA/l  AI’LILS  10  A  TKCC -PRE  E -f  Rf  E  -  ARf  F  PLAIE. 

C 

C  THE  VI3COELA31IC  MATERIAL  13  5M  ISD-112. 

C 

c  HIE  FoiLoiiiiiii  conririiiirG  arc  orriiiEO: 

C  El  -  YUdllG'S  liojlll  US  or  nA3E  KLAIE  (PSl) 

C  E2  =  'IOIIIIG'5  IIUDIIitlS  OF  VISCOFLASIIC  lArER  (PSl) 

C  Ei  =  VOUIIG’S  MOIHILUS  OE  COII3 1  RA  I  111  MG  LAYER  (I'Sl) 

C  G2  -  SMEAR  IU1LIULU3  OF  VISCUEIASIIC  IIAIEKIAL  (I'Sl) 

C  Mill  =  I'UISSOII'S  RATIO  OF  DASE  PIAIE 

C  1102  --  I'UISSOII'S  RAIIO  OF  VISCOEIASTIC  MATERIAL 

C  Ml  =  1M1CKMC3S  OF  BASE  PlAlE  (HI) 

C  II2  =  HIKIIIESS  or  VISCOELASTIC  I  AYER  (HD 

C  Mi  =  IMICKMLSS  OF  COII3  I  R  A  I  M 1  MG  (AYER  (HD 

C  MIOT  =  lOIAL  PLAIE  TMICKMESS  (HD 

C  KMOl  :  IIEM5ITY  OE  BASE  PIAIE  (  L  I)  E  -  5  EC  »  *  2-' 1  M  *  »  A) 

C  RM02  =  IILMSIIY  or  V  1  SCOE  L  A  S  )  J  C  MATERIAL 

C  KMUi  =  OEMSIIY  or  COIISTRAIMIMG  LAYER 

C  T  0  ,  FROM.  MROM,  M,ML  ,  ETFROL  .  SL  ,  511,  FROI  ,  8  C  ARE  COETEICIEMIS 

C  FOR  HIE  REUUCCn  FREQUEIICY  MOMOGRAM  EOUAllOUS 

C  FP  =  MOUAL  EREQUEMCY  OF  THE  UIIDAMPED  BASF  PLATE  (MZ) 

C  TCP  =  MODAL  FREQUEIICY  OF  HIE  COMPOSIIE  PLAIE  (HZ) 

C  ETA2  =  LOSS  FACTOR  OF  VlSCOELASllC  MAIERIAL 

C  EIAS  =  SYSTEM  LOSS  FACTOR 

C  KQR  =  HAVE  IIUMBFR  OF  PLATE 

C  GC  =  GRAVITAT lOMAL  COHSTAHT  (ni/SCC**2) 

C  T  5  TEMPERATURE  OF  V I  SCOFiL  AS!  1C  MATERIAL  (DEC  F) 

C  UP  =  MODAL  FREQUEIICY  OF  (UIOAMPED  PLATE  (RAD/SEC) 

C  FP  =  MODAL  FREQUEIICY  OF  UNDAMPED  PLAIE  (HZ) 

C  HCP  =  MODAL  FREQUEIICY  OF  DAMPED  PIAIE  (RAD/'SEC) 

C  FCP  =  MODAL  FREQUEIICY  OF  DAMPED  PLAIE  (HZ) 

C 

C  THE  UNITS  USED  IN  THIS  PROGRAM  ARE  LB,  INCH,  SEC,  DEGREES  F 
C 

C  UNDAMPED  MODAL  FREQUENCIES  ARE  FROM  A  FINITE  EIEMENT  ANALYSIS 
C  OF  HIE  FREE-FREE-FREE-FREE  PLAIE. 

C 

C  RESULTS  ARE  OUTPUT  TO  DATA  FILE  "lOSFCTR  DATA" 

C  PRIOR  TO  RUMNIIIG  THE  PROGRAM  TYPE  THE  COMMAIlD 
C  "FILEDEF  LOSFCTR  DISK  LOSFCTR  DATA" 

C 

C  FINITE  ElEMENI  MODAL  (REQUENCIES  ARE  INPUT  FROM  FllF. 

C  "PLirRiJ  DATA" 

C 

C  XXXXXXXXXXXXXXXKXXMXXXXXKXXXXXXXXXXKKXXXXXXXXXXXXXXXXXXKXXX 

C 

REAL  El  ,  C2,  Ei,  MUl  ,  IIU2,  111  ,  N2, 115,  RMOl  ,  RM02  ,  RHO  J,  III  01 

REAL  10, FROM, MROM, N, ML , ETFROL, SI, SN, FROL, C 

REAL  FP, FCP, ETA2, El AS,KQR,P1 ,GC, T 

REAL  T  R10,C1A210,SUIH  , SUU2 , SUB 5 , M 1 0 

REAL  1121  ,1151  ,G.C1  ,Al  PMRE.ALPNUI.  TIRE,  DIM,  DFIRE 

REAL  DEL  III,  EI(CUDC,SUDA,  SUD3,MP,IICP,  DENS,  SUD6  ,  SUB 7  ,  SUBS 

HHEGER  V 

DIMENSION  1IP(8) 

PI  =  S  .  XAT  AIK  1  .  ) 

C 

0Prtl(UNIT  =  10,Fn  E  =  'inSrCTR',STATU5='ni  D*  ) 

Ol’EIKUIII  m  1  1  ,  I  1  L  L  =  'PL  I  I  R'J'  ,  SI  A1U5  =  '01  O'  ) 

C 

C  ASSIGN  MMLRIAl  COMSIAIIIS  FOR  BASE  P I  A  T  F  ,  COIIS  I R  A  HI  I  MG  LAVER 
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C  MID  VIoCOriAriTIC 

c 

f  1  =  1 , 0  r  7 
1 1 1  =  U  .  7  "j  U 
HMOl  -  .  U'/AS 
iiiii  = .  ’ : 

RHU2= .035 

e:,m  .or? 

KllOi-  .  0968 
(^i:  =  386  . 

I1U2=.5 

c 

C  DCriNE  C0MSTAIIT3  FOR  NOMOGRAM  CQUAIIOII5 
C 

T0  =  1 OA  .  0 
rROM=2 . or9 
MP,OM  =  688  .96 
li-  .  275 
ML  =3 . 7 
ETFROL  =  1 . 08 
5L  =0  .  65 
5H^- . 55 
FROL  =5000 . 0 
C  =  2.5 
C 

C  DA3E  PLAIE  THICKNESS  LOOP 
C 

HI  =  0  .  37  5 
no  100  J  =  1 , 6 

IIRIIEL 10,703)  'BASE  PLATE  THICKNESS,  H1=',IU 
703  FORMAT! A27 , F6 . 3) 

HRI 7  E( 10, X) 

C 

C  VISCOELASTIC  THICKNESS  LOOP 
C 

H2=0.015 
DO  150  V=l,6 
C 

C  CONSTRAINING  LAYER  THICKNESS  LOOP 
C 

H3=0.0625 
DO  200  L=l,6 

IIRITE(10,702)  'VISCOELASTIC  LAYER  THICKNESS,  N2=',N2 
IIRI  lEC  10,») 

HRI rE(l0,702)  'CONSTRAINING  LAYER  THICKNESS,  H3=’,H3 
702  rORMAT(A35,F6 .6) 

HRITEC 10,X) 

HRI TEC  1  0,700)  'TEMP ' , 'MODE* . ' FCP ' ,  ’  E T A2  '  ,  ' ET AS ’ .  ' G2 ' 

7  00  FORMAT  CAS, 3X, A5, 6A1 5) 

C 

C  PLATE  MODE  LOOP,  READ  MODAL  FREQUENCY  AND  COMPUTE  HAVE  NUMBER 
C 

MTOT  =  HHH2>Hi 
DO  300  1=1,8 
READ! 1 1 , K )  UPC  I  ) 

C 

C  CALCULATE  MODAL  FREQUENCY  AND  HAVE  NUMBER  OF  BASE  PLAIE 
C 

SUB7  =  30R  I C  (  El  *CNT0rxx3)xGC)/'(  1  2  .  x(  1  . -NUl  KX2)  KRHDl  *111  )  ) 

KQR  =  1TPC  1  )/SUD7 
rp=HP( 1 )/(2 .*ri ) 
c 

C  TEMPERATURE  LOOP 
C 

T=30 , 0 

DO  600  K=1 , 15 
C 

C  CALCULATE  PROPERTIES  OF  VISCOELASTIC  FOR  GIVEN  TEMPERATURE  AND 
C  MODE 
C 

501  FR10=IOG10CFP)-C12.X(T-T0))/(525.<I-T0) 

F  R  =  1  0  .  C  (  R  I  0  ) 
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A  -  ( I  f;  1 0  - 1  in';!  0(  I  ROL )  )-'C 

sum  « ( (  SL )  sii  j t  ( SI  -sii i .  -sqr i  ( i .  i  Ax)r? . )  >  >/? . 

Cl  A21  U  =  l  (Hil  0(  C  I  ruol  )  I  SUUI 

c 

C  VISCCCLASIIC  mss  lAClUR 

c 

c  I A  2  m  0 .  • » ( r  I A  m  0 ) 

SU1)2  =  ?  .  »L(Uil  (U  llRUM.'tlL  ) 
suiis^i  .1(1  Rnn/rR)H'<ii 
111  0  =  L  OGl  0(  lie  )  «  SUD2/SUn5 

c 

C  VI5COELAS11C  SHEAR  MODULUS 
C 

G2m  0  .  x*  (Ml  0  ) 

C 

C  VISCOELASTIC  YOUUG'S  MODULUS 

r- 

E2  =  G2»2  .  XC 1  . ; HU2) 

C 

C  CALCUIATIOMS  FOR  DAMPED  PLATE  AT  TEMPERATURE  T  USING  R05S- 
C  KCRIII  ll-UHCAR  E;3UATIUIIS 
C 

G=G2/(E5*H5*H2XK0R) 

1121  =(H1  tll2)/2. 

1131  =112*  (  111  HI3  )/2  . 

Cl  =  C1  *111  x(  1  .  iG)  *G*E5xll5 
D  =  G''£1  A2*(  FI  *111  4  E  JXH5) 

AIPIIRE  =  G*E1XII1»C3*HJ)1(II31XM2.)><(C1*D*ETA2) 

ALPMIM  =  G*E1*IU  *E3*Hi>!(Mil*K2.  )kC1  *EIA2 
DRE  =  E1XII1XE2*II2*M31K(C1'*D)IEIA2) 

BIM  =  E1*II1  *E2*II2*HJI!<(C1XETA2-D) 
suns  =  2  .  XG* E2*M2*E 3*113X1121  XH5 1  *C 
DCLRE=SUB8*((1.-(ETA2*k2.))+(D*(2.*ETA2))) 
DrLIM=SUB8*(Cl*(2.*ElA2)-(D*(l.-(ETA2**2.)))) 

SUDA  =  (  12  ./(Cl  **2.  D**2  .  )  )*(  ALPHRE-BRE-DELRC) 

EMCUDE  =  (E1*(M1**3.  )  H  E3*(  113**3  .  HSUB-^ 

C 

C  MODAL  FREQUENCY  OF  DAMPED  PLATE 
C 

DF.NS  =  RH01*imRtl02*H2«H5»RH0.3 
SUD5=(EHCUDE*GC)/( 12.*(1 .-NU1**2. )*DENS) 

IICP  =  KqR*SQRl  (  SUDS) 

FCP  =  llCP/(2.  *ri ) 

C 

C  COMPARISON  OF  FP  AND  FCP 
C 

IF  (ADS( 1  .-FP/FCP) .LE.  0.1)  THEN 
GOTO  500 
ELSE 

FP=FCP 
GOIO  SOI 
ENDIF 
C 

C  COMPUTE  SYSTEM  LOSS  FACTOR 
C 

500  5UD6  =  (12./(C1**2.+T)**2.  )  )  M  (  Al  PH  I  IT- D I M- DE  L  I M ) 

CTAS  =  (  1  ./EIICUnC)*SUD6 
C 

C  PRINT  RESULTS 
C 

MRITEd 0,701 )  T, I, FCP, El A2. E T AS. G2 
7  01  rORMAT(5X,F7.5,2X,  I  2 , 3X,  <tE  1  5  .  A  ) 

C 

C  NEXT  temperature 
C 

T=T+5.n 

AOO  CONTINUE 
C 

C  NEXT  MODE 
C 

HR  1  TEC  10, *) 

300  CONTINUE 
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RCIIlMIKUlU  r  =  ll  ) 


c 

c  iiLxr  cniisj RAi 111 iitj  lArtR  iiiicKiir-'j-i 
c 

n.’i  =  ii2;i  n .  0625 

200  cdiiiiiiuf: 

c 

C  lltxr  VISCOELASIIC  LAYER  IH1CKIIL35 
C 

II2-H?'  n  .015 
150  COIIIll.'UC 
C 

C  NEXT  BASE  ELATE  TIIICKIIE5S 
C 

111  =HH  n  .  1  25 
100  cniiriiiuE 

Cl  osE(t)iiir  =  io) 

CLOGECUIlII^ll  ) 

EIIU 
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APPENDIX  B 


FORTRAN  PROGRAM  USED  TO  COMPUTE  MODAL  LOSS  FACTORS 
FOR  THE  DOUBLE  DAMPING  LAYER  CONFIGURATION  DESIGN 

This  program  uses  the  Ross— Kerwin-Ungar  equations  of  Chapter  3 
to  compute  estimated  modal  frequencies  and  loss  factors  for  the  double 
layer  configuration.  Modal  frequencies  for  an  undamped  plate  are  read 
from  a  data  file  and  estimated  modal  frequencies  and  loss  factors  for 
various  layer  thickness  combinations  in  the  double  layer  configuration  are 
output  to  another  data  file.  Material  properties  of  ISD-112  are  computed 
using  University  of  Dayton  data  and  curve-fitting  equations  to  the  reduced 
frequency  nomogram  [Ref.  9,12].  The  units  used  in  this  program  are 
pounds,  inches,  and  seconds. 


100 


on  n  n  000000000000000000000000000000000000000000000000000  r/r.^rj  noon 


rf;ut;i;An  iiiulyr 


XXXXXXKXKXHXXXMXXXMXXXyKXKXMtOIXKKKKIINKNXKIIXNXKKXMXKMItllXXXXKXXXKXtIIAXX 

nUG  t'f'OORAtl  CUMPUIER  TtlE  SYSll.ll  1050  TAriOR  AIID  MODAl  FREQUfriCY  01 
A  lUIOlllF  L'OIIO  I  liA  I  HI.  i)  LAYfR  V  I  OC:  U  E  L  AS  I  I  C  A  L  L  Y  DAlirEl)  I'LAIE.  IHE  LOOS 

fAcioir,  Al'c  rniiruici)  run  a  orrcirjc  ri a  1  r./ haih' im;  layer 

com  J  I.URA  I  1  oil  Alll)  OVER  A  T  Elll' L R A  1  URL  KAlir.E  01'  iO  lOO  UEGREES 
r  AHLimE  I  r . 

Mils  EROORAH  APPLIES  TO  A  E  RE  P  -  E  R  EE -E  R  f  F  -  f  R  EE  I'LAIE  ADI)  MIC  UDOAItPEU 
IlnpAL  IKCOULNCir.S  ARE  UCICRnilM;l)  f  ROII  Ml/1  IE  ELEHEM1  AIIALYS15  Alll) 

ARE  READ  ll/TO  llllS  PROGRAM  PROM  PILE  'I’LlIRi)'. 

THE  VISCOEIASTIC  MATERIAL  IS  3M  lSD-112.  VISCOELASTIC  MATERIAL 
DATA  IS  PROM  UMIVERSIIY  OF  DAYIOM  RESEARCH  IIISIIIUIE 

THE  UIIIIS  USED  III  THIS  PROGRAM  ARE  LD,  INCH,  SEC,  AIID  DEGREES  P. 

LOSS  PACTORS,  DAMPED  PLATE  MODAL  FREOUEHCIES,  AIID  ISD-112  PROPERTIES 
ARE  COMPUIEU  III  SUIIROUTIHE  ’RKU*  f(IR  EACH  ll-IH  COIIS I  RA I II E  D  LAYER 
SYSTEM.  UPUri  COMPLEIIOII  OF  EACH  lEMPCRAIlIRE  COMPUIATIUH  THE  SYSItM 
LOSS  FACTOR  AIID  C OR R E SPUI/D 1 IIG  DAMPED  PLAIC  EREQUEIICY  ARE  tlRlTTEII  10 
r  1  L  E  '  T IIOLYR  DAT  A  '  . 

THE  rni  LOMitiG  cocrr ici EiiTs  are  oErinED: 

El  =  YOlIHG'S  MODULUS  OF  HASP  PLAIE  (PSI) 

Pi  =  EQUIVAlElir  YOUIIG'S  MODUl  US  OP  N-IH  COIISIRAIIIRD  LAYER  SYSUM 
EiPRM  =  YOUIIG'S  MODULUS  OP  THE  COIIS  1  RA  1 II 1  IIG  LAYER  IH  M-Tll  LAYER 
EHCUDE  =  EOUIVALEIIT  STIFFNESS  OF  COIIS  I  RA  I  MED  LAYER  5YSIEM 
ETAS  =  SYSTEM  LOSS  FACTOR 

EFA3  =  LOSS  FACTOR  OF  M-TM  COHSTRAIMED  LAYER  SYSTEM 

ETA3PM  =  LOSS  FACTOR  OF  COIIS  I  RA  1  MI  IIG  LAYER  Ill  HIE  ll-IM  LAYER  (  =  0 

TCP  =  FREDUEHCY  OF  IHE  DAMPED  PLATE  (HZ) 

FP  =  FREQUEIICY  OF  UNDAMPED  PLATE  (112) 

GC  =  GRAVITATIONAL  CONSTANT 

111  =  THICKNESS  OF  DASE  PLAIE 

H2  =  THICKNESS  OF  1ST  VISCOELASTIC  LAYER 

113  =  THICKNESS  Of  N-TH  COIIS  I  RA  I NE  D  LAYER  SYSTEM 

liA  =  TOTAL  DAMPED  PLATE  TMICKHESS 

tllPRM  =  THICKNESS  OF  BASE  LAYER  IN  H-IM  LAYER 

H2PRM  =  THICKNESS  OF  VEIl  III  H-TH  LAYER 

H3PRM  =  TMICKHESS  OF  COHSTRAIHItlG  LAYER  IN  M-TH  LAYER 

KQR  =  HAVE  HUMBER  OF  UNDAMPED  PI  ATE 

IIUI  =  POISSON'S  RATIO  OF  BASE  PLAIE  AMD  CONST  RAINING  LA'.LRS 

11112  =  POISSON'S  RATIO  OF  VISCOELASTIC  MATERIAL 

NU2PRM  =  POISSON'S  RATIO  OF  VEM  IN  ll-TII  LAYER 

RMOl  =  DENSITY  OF  BASE  PLAIE  AND  COIIS  I  KAl  N1 NG  LAYERS 

RN02  =  DENSIIY  OF  VISCOELASIIC  MATERIAL 

RN02PM  =  DENSITY  OF  VEM  IN  ll-TH  LAYER 

nit03  =  DENSITY  OF  N-TH  LAYER  (=  RHOl  ) 

RN03PM  =  DENSITY  OF  N-TN  LAYER  CONSTRAINING  LAYER 
T  =  TEMPERATURE  VARIABLE 

HP  =  FREQUENCY  OF  UNDAMPED  PLAIE  (RAD/SEC) 

C,ETFROl  ,FROL,FRnM,ML,MROI1,N.SN,SI  .TO  =  COEFnCIENIS  FOR  THE 
REDUCED  FREQUENCY  NOMOGRAM  EQUAII0N5 

xxxxxxxxxxkkxxxkxxxxxxxxxxxxxnkkxxxxkxxxxxxxxxxxxxxxxxxxkxxxxxxxx/cxxxx 

REAL  C, El , E3, E3PRM, EHCUDE, El  AS, ET A3, ET A 3P11 , E T P ROL , TCP, PP, PROl 
REAL  PROM,  GC.  Ill  ,  112 , N 3 , IIA  ,  H 1  PRIl ,  II2PRM ,  N3PRn ,  KQR  , Ml  ,MRUM,  ll.NUl 
REAL  NU2,NU2Pn,  PI  ,  RNO 1 ,  RII02 ,  RII02PM,  SN.  SL  ,  SUDl  ,  I  .  TO.IIP 
REAL  RN03,RN03PM 

DIMENSION  IIP(  17  ) 

PI  =0  .  oxaIAIK  1 .0) 

OPEIK  IINIT  =  10,  PII  E  =  'PI  TPRQ'  ,  S I A 1115=  '  01  D '  ) 

0PEII(UIIII  =  11,PILE-'  IMULYR',SIAIUS='OLD*  ) 

ASSIGN  LAYER  THICKNESSES 
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Ill  =  0  .  -i 
1 1  r  -  u  0 5 


M  1  1 

i;ii 

0 

.  0  T  3 ; 

5 

1121 

on 

0 

.015 

1131 

i.'ii 

0 

09  37 

5 

113  = 

III! 

■p 

1 1  1121' 

Kin  II3PRI1 

1 1  'I  = 

Ill 

II 

2  1  113 

c 

C  A5Sir.ll  llAiriilAL  Clllir.  I  All!  5 

c 

C 1 ^ I . 0C7 
KIIUI  =  0 . 0968 
mil  =0  ,  L5 
RII02  =  0  .  0i5 

mi2  =  0 . 99 

RII02f'ri  =  0 , 035 
111120,1:0  .  5 
RII03:0 . 0968 
E3f'RM=l  .0C7 
RllOSrn^O  .0968 
El  A3l'n^0  .  0 
ri A3  :  0.0 

c 

GC:386  .  0 
C 

C  DEEIIIC  VISCOELASTIC  COIISTAIfTS  FOR  tlOllOGRAM  EQUAllOtlS 
C 

T  0  =  I  0  A  .  0 
rRon=2 . OCA 
|;R()I1:688  .99 
11  =  0.275 
IIL  =8  .  7 
ETFROL  =  1.08 
SI.  =  0.95 
5ll  =  -0 . 55 
FROL  =5000 . 0 
C  =  2.5 
C 

C  IIRITE  PLATE  CHARACTERISTICS 
C 

ITRITC(  1  1  ,  700  )  '111  ='.111 
HRITEdlWOO)  M12  =',H2 
HRITEdlWOO)  'MIPRIME  =',H1PRM 
llRIIEdl,  7  00  )  •II2PRIME  =',II2I’RM 
IIRI  I  Ed  1  ,  700  )  'M3PRJI1E  =',M3rRll 
IIRI  IE(  Id*) 
mu  TE(  11  ,  X  ) 

700  rORnAT(A12,r6.5) 

IIRI  T  Ed  d  7  01  )  'I  EMP'  ,  'MODE'  ,  TCP*  ,  '  El  AS' 

701  FORMATCAS, 3X,A5,2A15) 

HRITEddx) 

C 

C  PLATE  MODE  LOOP 
C 

DO  100  1  =  d 1 7 
REAUCIO.X)  IIP(I) 

C 

C  CALCULATE  MODAL  FREQUEIICY  AMD  HAVE  IIUMDER  OF  UllUAMrED  PIAIE 
C 

SUDdSORTdElxdlAx.xSlxCO/dZ.xd  .-IIU1kk2)*RII01*II9)) 

KOIdlll’d  )/SUDl 

ri’-lll’(  I )/(  2  .  XPI  ) 

C 

C  TEMPERATURE  LOOP 
C 

T=30 . 0 

DO  200  K=d  15 
C 

C  COMPUTE  SYSTEM  LOSS  FACTOR  AMD  SIIFFMCSS  TOR  M-IH  LAYER 
C 

CALL  RKUdllPRM,  HCPRi:,  HiPR:!,  tdRtlOd  MUl  ,  in;2.’il,  RII02II1,  E3PRIi,  RllidPli, 
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nonononnonoonnoorononnonnnoooonoo  nnno  o>-  nnnrj  nnn  r,  on 


Cl  I  A  ’.III,  I  ,  I  I' ,  KIJK,  Cl  AS.  tllCIJUL,  rcr,  I  0  ,  IKOIl,  UKOd,  II,  m  ,  cn  KOL  ,  SL  ,  Sll. 

Cl  icm  ,  c  ) 
c 

C  fOlIVlKI  RCStiLTS  CROn  ll-lll  lAYLR  CAlCUl  A I  I  011  10  lUIAL  I'lAlE 

c 

r I a:=  n  .  u 

I  R  =  lire  I  )  /  ( .  KPn 
c  ir^rcr 

r  i  - 1  lu  einc-'  ( ii:«  *3 ) 

c 

c  cunruiE  SYSicii  loss  eacior  aiid  ercqueiicy  ruK  ioiai  riaie 

c 

CALI  RKU(  ll  1  ,  112,  M5,  Cl  ,  KIlOl  .  Mill  ,  MU2,  RII02  ,  E ,  R 110 ,  E  T  A  3  ,  T  ,  E  P  ,  K.Q  R  ,  E  1  A  S  , 
CCIiCUDE  ,  ECP  ,  I  0  ,  EROIl,  IIROM,  tl.IlL  ,  C  1  E ROL  ,  SC  ,  Sll ,  E KOL  ,  C  ) 

ERIIIT  RESULTS 

"HI  TEC  1  I , 702  )  T, I, ECP, ETAS 
EORMAT  (  5X,F7.3,2X;l2,5X,2E15.‘i) 

NEXT  TElirERAIURC 


T  =  T  I  5  .  n 

00  coiiiinuc 
NEXT  none 

IIRI  TEC  1  1  ,  *  ) 

00  COMllllUE 

CLOSECtmiT=ll) 

CLOSECUniT=10) 

EIID 

XXXXXXXXXXXKXK.'OCXXKXXXXKXXKXXXXXKKXXXKXXXXXXX'O'XXXXKXXKXXXXKXXK******* 

X.XXXXXXXXXXXXXXXXXXKXXXXXXXXXMMKXXK.XKKXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXX 

SUDROUTIME  RKUC  HI ,  H2  ,  M3,  El ,  RHOl ,  tIUl ,  IIU2,  RH02 ,  E5  ,  RIIO  5  ,  E  T  A  3  .  T  ,  FP  , 
CKQR,ETAS,EIICUBE,rCP,TO,FROM,MROM,fl,ML,ETrROL,SL,SH,FROL  ,C) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

THIS  SUBROUTINE  CALCULATES  VISCOELASTIC  PROPERTIES  BASED  ON  THE 
UNIVERSITY  OF  DAYTON  REDUCED  FREQUENCY  NOMOGRAM  EQUATIONS,  AND  TIICN 
CALCULATES  PLATE  STIFENESSES  AND  LOSS  FACTORS  BASED  ON  THE  ROSS- 
KERHIN-UNGAR  EQUATIONS. 

THE  FOLLOIIING  ADDITIONAL  VARIABLES  ARE  DEFINED  FOR  USE  IN  THIS 
SUBROUTINE- 

A  =  COEFFICIENT  FOR  NOMOGRAM  EQUATIONS 

ALPHIM, ALPHRE  =  IMAGINARY  AND  REAL  COMPONENTS  OF  COEFFICIENT 
ALPHA  IN  THE  RKU  EQUATIONS 

BIM.BRE  =  IMAGINARY  AND  REAL  COMPONENTS  OF  COEFFICIENT  'B'  IN 
THE  RKU  EQUATIONS 

Cl.D  =  COEEFICIENIS  FOR  RKU  EQUATIONS 

DEL1M,DELRE  =  IMAGINARY  AND  REAL  COMPONENTS  OF  COEFFICIENT 
DELTA  IN  THE  RKU  EQUAIIONS 

DENS  =  COMBINAIION  OF  MATERIAL  DENSITIES  USED  TO  COMPUTE  THE 
FREQUENCY  OF  THE  DAMPED  PLATE  AND  DAMPING  LAYERS 
EllCUBE  =  EQUIVALF.NI  STIFFNESS  OF  DAMPED  PLATE  AS  COMPUTED 
USING  RKU  EQUAIIONS 

ETA2  =  LOSS  FACTOR  OF  VEM  COMPUTED  IN  NOMOGRAM  EQUATIONS 
EIA210  =  L0GI0(ETA2) 

EIAS  =  SYSTEM  LOSS  FACTOR  COMPUTED  BY  RKU  EQUATIONS 
FCP  =  MODAL  FREQUENCY  OF  DAMPED  PLATE  < N2 ) 

FP  =  MODAL  FREQUENCY  OF  UNDAMPED  PLAIE  (HZ) 

ER  REDUCED  FREQUENCY  OF  VEM 
FRIO  =  LOGin(FR) 

G2  =  SHEAR  MflDUmS  OF  VEM 

1121,1131  =  RKU  EQUAIIUN  COEEF  1 C  1  EN T S 

MIO  =  LI)G1U(G2?  AS  COMPUIED  BY  NUMUGRAM  EQUAIIONS 
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c  iici’  =  rRi.gunicr  or  DAnrco  I'LAirr  (rad/gco 

c 

REAl  A,AI  n!int,AirMRC,nil1,BRC,C,Cl,D,[)LLni,UEI  RE,I)rnS,ri,E2,Fi 
REAL  Lliconc,  E  r  A?,  L  I  A5,  El  A?1  0,  F  I  A5,  E  II  ROL  ,  rcr,  rr,  ER,  TRIO,  EROL  ,  (  ROM 
RCAl  G.r,2,l^C,IU  ,112,  Hi,  1(21  ,  Hi  1  ,  KQR  ,  H 1  0  ,  HE  ,  HRGM  ,  II ,  HU  1  ,HU2,n  ,  RIUH 
REAl  R|l02,R'll05,jH.5l.SUIIl,5llD2,5UDi,SUD'<,5UD5,3Uli6,GUn7,SUD8 
REAL  GUb9,GU1110,5UUll,r.IO,llCP 
C 

n -A . OXAIAIK  1.0) 

GC=i86 .0 
C 

C  CALCULATE  EROrERTIES  OF  VEH  FOR  GIVEH  TEMPERATURE  AMD  MODE 
C 

501  rR10=LOG10(rP)-(12.x(T-TO))/(525.^T-T0) 

FR=10.xx(rR10) 

A-( f  R1 O-lOGI 0( FROL  )  )/C 

SUDl=CXtCSL«SH)xA«(GL-3H)X(l.-SQRI(l.«AXx2)))/2. 

ErA210  =  LOG10(Err  ROD  *5001 
C 

C  VISCOELASTIC  LOSS  FACTOR 
C 

CTA2=10.xx(ETA210) 

C 

C  VISCOELASTIC  SHEAR  MODULUS 
C 

SUD2  =  2  .  XLOGlOCMROn/IIL  ) 

SUBi=l .+(FEOM/FR)xxH 
M10=LOG10(HL )+SUD2/SUD3 
02  =  10  xxdllO) 

C 

C  VISCOELASTIC  YOUMG'S  MODULUS 
C 

E2=G2X2 . »( 1 , +MU2) 

C 

C  CALCULATIONS  FOR  DAMPED  PLATE  USING  RKU  EQUATIONS 
C 

G=G2/(E3*H3XH2KKQR) 

H21  =  <HHH2)/2. 

H31=H2+(Hl+H3)/2. 

Cl=ElxlUx(  1  .+0)  +  G*E3)*H3!<(l.-ETA2»<£TA3) 

D  =  G!<E1*N1  xETA2  +  CxE3)*M3*(  ETA2»ETA3) 

5UDA  =  GxE1kM1xE3!<H3x(H31xx2) 

ALPHRE^SUDAXCCl x( 1 . - E T A2x ET A3 ) + Dx ( ET A2+ E T A3 ) ) 
ALPHIM=SUBAx(C1X(ETA2tETA3)} 

SUB5=E1XN1XE2XN2XH31 

BRE=SUD5X(C1+DXETA2) 

BIH  =  SUD5X(C1XF.TA2-D) 

SUB6=2. XGXE2XN2XE3XH3XH21XH31 
SUD7=1 .-2.xETA2xETA3-(ETA2KX2) 

SUn8  =  2.’<ETA2  +  ETA3-(ETA2xx2)xErA3 
DELRE  =  SUD6x(ClxSUB7  tDxSUBB ) 

DEL1M=SUB6X(CIXSUD8-0XSUD7 ) 

C 

SUBO^I 12./(C1 xx2( DXX2) )X( ALPHRE-BRE-DELRE) 

ENCUDE=Elx(|llxx3) I E3X(H3xx3)tSUB9 
C 

C  MODAL  FREQUENCY  OF  DAMPED  PLATE 
C 

DENS  =  RN01xM1*R1I02xN2  +  RH03xM3 

SUB10  =  L  EIICUDExGC)/(  1 2  .  x(  NUI  Xx2  )  xDENS  ) 

HCP=KqRXSQRT(SUD10) 

FCP  =  MCP/'(2.XPI) 

C 

C  COMPARISON  OF  FP  AND  FCP 
C 

IF(  AD5(  1  .-FP/FCP)  .LE.  O.IO)  TMr;N 
GOIO  500 

ELSE 

FP=FCP 
Goro  501 
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Lfini  F 


c 

C  COMI'UIE  SYGIEtl  IUS5  FACIOR 
C 

50  0  5UB1  1  -  (  1  2  .  /  (  Cl  X  I'.ix?)  )*(  Al  PHiri  DIM  -  DEL  IM) 

EIAS=(L./EHLUnE)x((Ei»li5xElAiItSUDll) 

C 

CfID 


105 


APPENDIX  C 


DESIGN  DRAWINGS  FOR  THE  MACHINING 
OF  THE  FLOATING  ELEMENT  AND 
POCKET  PLATE  CONFIGURATIONS 

The  drawings  shown  in  Figures  C.l  and  C.2  were  used  to  machine 
the  pocket  plate  and  floating  element  plate  used  in  the  experiments.  These 
drawings  are  included  to  show  the  relation  of  the  pocket  for  the  ISD— 112  to 
the  cover  plate  and  how  the  viscoelastic  was  protected  from  the  heat  of 
welding. 
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Figure  C.l.  Design  Drawing  of  the  Pocket  Plate  Configuration. 
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Figure  C.2.  Design  Drawing  of  the  Floating  Element  Configuration. 
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APPENDIX  D 


REPRESENTATIVE  MSC/NASTRAN  DATA  DECK  FOR 
THE  DAMPING  CONFIGURATIONS 

This  data  deck  was  used  to  compute  the  modal  frequency  response  of 
the  single  layer  damping  configuration  and  is  a  representative  sample  of 
the  NASTRAN  decks  used  for  the  other  finite  element  models.  The  values 
in  the  damping  table  are  from  a  cur\'e  fit  to  the  modal  loss  factors  estimated 
from  the  modal  strain  energ>’  method.  Since  the  data  deck  for  the  normal 
mode  and  modal  strain  energy  extraction  is  \drtually  identical  to  this  decK, 
the  Case  Control  deck  commands  for  the  normal  mode  analysis  are 
included,  but  are  commented  out. 

The  OUTPUT  request  provides  data  for  an  x-y  plot  of  the  modal 
frequency  response. 

The  units  used  in  this  deck  are  pounds,  inches,  and  seconds. 
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Id  Slagle, ifr 
sol  30 
$ 

$  THIS  DECI  IS  TO  COMPUTE  TEE  MODAL  FIEQUEICT  lESPOISE  OF  A  15x45 

$  IICE  ALUMIIUM  PLATE  HITE  A  COISTIAIIED  VISCOEUSTIC  DAMPIIG  LAYEI. 

$  TEE  MODEL  EAS  252  ELEMEITS  HITE  64  ELEMEITS  II  EACE  LATEI.  THE  QUAD 
$  ELEMEITS  AIE  OFFSET  FEOM  TEE  lEl  ELEMEITS  AS  SUGGESTED  BY  TEE 
$  LITEIATUIE. 

$ 

$  THE  MATEIIAL  PIOPEITIES  OF  ISD-I12  AIE  FIOM  THE  3M  COlPOBATIOI. 

$ 

$  TEE  NODAL  LOSS  FACTOIS  II  TEE  DAMPIIG  TABLE  ABE  FROM  A  CUIVEFIT 

$  TO  TEE  SET  OF  MODAL  LOSS  FACTOIS  COMPUTED  FIOM  THE  MODAL  STIAII  EIEIGY 
$  METHOD. 

$ 

$  TEE  UIITS  USED  II  THIS  DECI  AlE  POUIDS.IICHES,  AID  SECOIDS 

tlie  3000 
ceod 

title  :  MODAL  FIEQUEICY  lESPOISE  /  200  HZ  /  3M 

■ethod  =  1 

spc  :  1 

dload  :  10 

Irequeocy  :  10 

sdaiplag  =  101 

set  111  :  95 

svector  *  all 

acceleratloD(plot, phase)  ^  111 

output(iyplot) 

xyprlDt  acce  /  95(t3) 

$ 

J  ttUIttlttiliit 

f  TEE  FOLLOHIIG  LUES  AIE  THE  CASE  COITIOL  DECI  CARDS  FOR  TEE 
$  lOIMAL  MODE  EITlACTIOl  AID  STIAII  EIEIGY  lEQUEST 
$  letbod  :  1 
I  spc  =  1 
$  set  10  :  all 
$  set  11  =  169, thru, 252 
I  ese  :  11 
$ 

$  . . . 

BEGII  BULI 

$  TITLE  =  SIIGLE  LAYEI  HITE  QUAD  OFFSET 
$  DATA  DEC!  PRODUCED  BY  PATIAS  VERSIOI  2.0:  24-I0V-89  06:31:46 


GRID 

1 

45.0000  15.0000  0.51500 

GRID 

2 

45.0000  15.0000  0.50000 

GUO 

3 

45.0000  12.5000  0.51500 

GRID 

4 

45.0000  12.5000  0.50000 

GRID 

5 

45.0000  10.0000  0.51500 

GRID 

6 

45.0000  10.0000  0.50000 
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GIID 

7 

45.0000 

7.50000 

0.51500 

GIID 

6 

45.0000 

7.50000 

0.50000 

GIIO 

9 

45.0000 

5.00000 

0.51500 

GIID 

10 

45.0000 

5.00000 

0.50000 

GIID 

11 

45.0000 

2.50000 

0.51500 

GIID 

12 

45.0000 

2.50000 

0.50000 

GIID 

13 

45.0000 

0. 

0.51500 

GIID 

14 

45.0000 

0. 

0.50000 

GIID 

IS 

41.7857 

15.0000 

0.51500 

GIID 

16 

41.7857 

15.0000 

0.50000 

GIID 

17 

41.7857 

12.5000 

0.51500 

GIID 

18 

41.7857 

12.5000 

0.50000 

GIID 

19 

41.7857 

10.0000 

0.51500 

GIID 

20 

41.7857 

10.0000 

0.50000 

GIID 

21 

41.7857 

7.50000 

0.51500 

GIIO 

22 

41.7857 

7.50000 

0.50000 

GIID 

23 

41.7857 

5.00000 

0.51500 

GIID 

24 

41.7857 

5.00000 

0.50000 

GIID 

25 

41.7857 

2.50000 

0.51500 

GIID 

26 

41.7857 

2.50000 

0.50000 

GIID 

27 

41.7857 

0. 

0.51500 

GIID 

28 

41.7857 

0. 

0.50000 

GIID 

29 

38.5714 

15.0000 

0.51500 

GIID 

30 

38.5714 

15.0000 

0.50000 

GIID 

31 

38.5714 

12.5000 

0.51500 

GIID 

32 

38.5714 

12.5000 

0.50000 

GIID 

33 

38.5714 

10.0000 

0.51500 

GIID 

34 

38.5714 

10.0000 

0.50000 

GIID 

35 

38.5714 

7.50000 

0.51500 

GIID 

36 

38.5714 

7.50000 

0.50000 

GIID 

37 

38.5714 

5.00000 

0.51500 

GIID 

38 

38.5714 

5.00000 

0.50000 

GIID 

39 

38.5714 

2.50000 

0.51500 

GIID 

40 

38.5714 

2.50000 

0.50000 

GIID 

41 

38.5714 

0. 

0.51500 

GIID 

42 

38.5714 

0. 

0.50000 

GIID 

43 

35.3571 

15.0000 

0.51500 

GIIO 

44 

35.3571 

15.0000 

0.50000 

GIID 

45 

35.3571 

12.5000 

0.51500 

GIID 

46 

35.3571 

12.5000 

0.50000 

GIID 

47 

35.3571 

10.0000 

0.51500 

GIIO 

48 

35.3571 

10.0000 

0.50000 

GIID 

49 

35.3571 

7.50000 

0.51500 

GIID 

50 

35.3571 

7.50000 

0.50000 

GIIO 

51 

35.3571 

5.00000 

0.51500 

GIID 

52 

35.3571 

5.00000 

0.50000 

GIID 

53 

35.3571 

2.50000 

0.51500 

GIIO 

54 

35.3571 

2.50000 

0.50000 

GIID 

55 

35.3571 

0. 

0.51500 

GIIO 

56 

35.3571 

0. 

0.50000 

GIID 

57 

32.1429 

15.0000 

0.51500 

1 1 1 


6IID 

58 

32.1429  15.0000  0.50000 

GIID 

59 

32.1429  12.5000  0.51500 

GIID 

60 

32.1429  12.5000  0.50000 

GIID 

61 

32.1429  10.0000  0.51500 

GIID 

62 

32.1429  10.0000  0.50000 

GIID 

63 

32.1429  7.50000  0.51500 

GIID 

64 

32.1429  7.50000  0.50000 

GIID 

65 

32.1429  5.00000  0.51500 

GIID 

66 

32.1429  5.00000  0.50000 

GIID 

67 

32.1429  2.50000  0.51500 

GIID 

68 

32.1429  2.50000  0.50000 

GIID 

69 

32.1429  0.  0.51500 

GIID 

70 

32.1429  0.  0.50000 

GIID 

71 

28.9286  15.0000  0.51500 

GIID 

72 

28.9286  15.0000  0.50000 

GIID 

73 

28.9286  12.5000  0.51500 

GIID 

74 

28.9286  12.5000  0.50000 

GIID 

75 

28.9286  10.0080  0.51500 

GIID 

76 

28.9286  10.0000  0.50000 

GIID 

77 

28.9286  7.50000  0.51500 

GIID 

78 

28.9286  7.50000  0.50000 

GIID 

79 

28.9286  5.00000  0.51500 

GIID 

80 

28.9286  5.00000  0.50000 

GIID 

81 

28.9286  2.50000  0.51500 

GIID 

82 

28.9286  2.50000  0.50000 

GIID 

83 

28.9286  0.  0.51500 

GIID 

84 

28.9286  0.  0.50000 

GIID 

85 

25.7143  15.0000  0.51500 

GIID 

86 

25.7143  15.0000  0.50000 

GIID 

87 

25.7143  12.5800  0.51500 

GIID 

88 

25.7143  12.5000  0.50000 

GIID 

89 

25.7143  10.0000  0.51500 

GIID 

90 

25.7143  10.0000  0.50000 

GIID 

91 

25.7143  7.50000  0.51500 

GIID 

92 

25.7143  7.50000  0.50000 

GIID 

93 

25.7143  5.00000  0.51500 

GIID 

94 

25.7143  5.00000  0.50000 
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